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ABSTRACT 


Experimental results of comparative screening studies of candidate molten 
carbonate salts as phase-change materials (PCM) for advanced solar-thermal 
energy storage applications at 540° to 870°C (1004° to 1600°F) and steam-Rankine 
electric generation at 400° to 540°C (752° to 1004°F) are presented. Alkali 
carbonates are attractive as latent-heat storage materials because of their 
relatively high storage capacity and thermal conductivity, low corrosivity, 
moderate cost, and safe and simple handling requirements. Salts were tested 
in 0,1 kWhr lab-scale modiiles and evaluated on the basis of discharge heat flux, 
solidification temperature range, thermal cycling stability, and compatibility 
with containment materials. The feasibility of using a distributed network of 
high-conductivity material to increase the heat flux through the layer of 
solidified salt was experimentally evaluated. The thermal performance of an 
8 kWhr thermal energy storage (TES) module containing LiKCOs remained very 
stable throughout 5650 hours and 130 charge/discharge cycles at 480° to 535°C 
(896° to 995°F) . A TES utilization concept of an electrical generation 
peaking subsystem composed of a multistage condensing steam turbine and a TES 
subsystem with a separate power conversion loop was defined. Conceptual de- 
signs for a 100 MWg TES peaking system providing steam at 316°, 427°, and 
454°C (600°, 800°, and 850°F) at 3.79 X 10*^ Pa (550 psia) were developed and 
evaluated. Areas requiring further investigation have also been identified. 


PBECEDINQ page bumk not FUJMEO 


INSTITUTE 


0 F 


GAS 


TECHNOLOGY 




11/79 


9100 


EXECUTIVE SUMMARY 

The objectives of this TES program were to a) define and identify a 
high-temperature TES concept that can effectively utilize the high-quality 
thermal energy in a steam-electric power generating facility, b) investigate 
and test candidate TES salts, primarily carbonates, which are applicable to 
energy storage in both the superheat region (temperature range of 343“ to 
39'^°C, 649" to 750°F) of conventlc<nal electric-utility, water-Rankine power 
plants and in TES systems for solar-thermal power system applications in 
the temperature range of 538“ to 871“C (1000“ to 1600“F) , and c) evaluate 
means of increasing the heat flux during thermal discharging. The work 
performed to achieve these objectives is described in two tasks. 

• Task 1 System Definition 

This section presents a TES peaking system concept to apply in fossil 
fuel-fired steam-electric plants. It discusses the design criteria used to 
develop the concept, the necessary operating characteristics from the TES 
unit, and the economics and general status of the technology. The principal 
features of the TES peaking system are 1) good system efficiency, 2) minimum 
degradation of the charging energy in the TES unit, and 3) minimum impact on 
the primary plant's operating procedures. 

The results of the system studies performed here indicate that a phase- 
change TES material for application in a steam-Rankine cycle should have a 
transformation temperature less than 350“C (662“F). 

• Task 2 TES Module Testing and Evaluation 

Three engineering-scale modules were constructed and operated under 
Energy Research and Development Administration (ERDA) Contract No. EY-76-C-02-2888 . 
The engineering-scale units were 0.305 m (1 ft) in outer diameter (OD) and 
0.457 m (1.5 ft) high with a 0.0508-m (2-in.) diameter single-pass heat 
exchanger (HX) tube. They contained LiKCOs salt, which has a melting point 
of 505°C (940°F) . One of these units, Unit No. 2, was cycled 30 times, 
which represented 1400 hours of operation. As a result of the already 
accumulated operating time, the unit was restarted and operation was continued 
under this contract to obtain system endurance data. This unit was then 
operated and discharged for a total of 5650 hours and 129 charge/discharge 
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cycles, with no significant losses of system performance. The LiKC 03 salt and 
AISl 316 stainless-steel containment and HX exhibited very good physical and 
chemical stability. 

Salts were also Investigated in this subtask for their applicability to 
latent heat storage in the low 343“ to 399“C (650° to 750°F) temperature range, 
and the 538° to 871°C (1000° to 1600°F) range for solar-thermal applications. 
Emphasis was given to salt compositions based on Li 2 C 03 , Na^COs, ^2^03* 
from which congruent melting carbonate mixtures are ayailtible covering the 
temperature range 397° to 898°C (747° to 1648°F). From the work completed in 
Task 1, however, it was indicated that phase-change TES materials in a 
steam-Ranklne cycle should also have transformation temperatures less than 
350°C (662°F) . Because this temperature is lower than the lowest melting 
carbonate eutectic known, and subsequent work was to be primarily with carbonates, 
the major experimental effort was at this point redirected to solar-thermal 
systems applications in the temperature range 538° to 871°C (1000° to 1600°F). 

These salts were tested in laboratory-scale units that were 0.0762 m (3 in.) 
in OD and 0.127 m (5 in.) high with a 0.0127-m (0.5-in.) diameter Internal 
HX tube. The thermal performance of the laboratory-scale units was 
determined during charge/discharge cycles, using high-velocity ambient air 
as the working fluid. From this study, pure Li 2 C 03 and Na 2 C 03 emerged as 
strong congruent melting candidates for high-temperature solar-thermal applications, 
with a (Li-Na-K) 2 G 03 eutectic for the lower, s team-Rankine applications. Two 
alkall/alkallne earth carbonate mixtures, Na 2 C 03 -BaC 03 and Li 2 C 03 -CaC 03 , also 
exhibited extremely stable performance in the 600° to 750°C (1112° to 1382°F) 
range. The thermal response of these lab-scale storage units agreed well with 
that predicted by the heat-transfer model developed at the Institute of Gas 
Technology (IGT) under the previous ERDA contract. 

At the flow rates investigated, the rate of discharge was limited by the 
transfer of heat between the salt and HX through the zone of solidified salt 
formed on the HX during discharge. Enhancement of this heat conduction was 
Investigated and accomplished by Installing a continuous, reticulated aluminum 
matrix into the (Li-Na-K) 2 CO 3 salt and by testing an off-eutectic, 81.3 wt % 

NazCOs — 18.7 wt % K 2 CO 3 incongruent melting slush. A porous carbon-based 
material tested as a heat conduction enhancement material was found to be 
imcompatible with a Na 2 C 03 -K 2 C 03 melt at 710°C (1318°F). 
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A heat transfer model, originally developed at XGT in 1976 under RRDA 
Contract No. EY-76-C-02-2888, was also used in the scale-up of a TES subsystem 
conceptual desxgn relevant to large, solar-thermal power system applications. 

A conceptual design based on use of Li 2 G 03 PCM for a 50 MW^h system utilizing 
air charge/discharge fluid was completed. 

It is recommended that further investigations focus on salt/containment 
compatibility using high-temperature alloys, development of cost-effective 
containment and insulation materials, long-term testing of larger units 
(>10 kWhr,-j^) utilizing actual working fluids, continued cost/benefit analyses 
of thermal conductivity enhancement (TCE) materials, and heat-transfer modeling 
studies. In addition, interactlou with ongoing TES systems design and 
analysis studies is necessary to provide effective feedback to salt/HX 
research and development efforts. 
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DETAILED DESCRIPTION OF TECHNICAL PROGRESS 
Task 1 . System Definition 
Introduction 

Energy demands from baseload electrical generating plants vary tliroughout 
any given day, producing a variable demand on generating equipment. Because of 
this variation, a portion of a utility’s baseload generation, which is generally 
coal- and nuclear-fueled, may be unused during low-demand periods. With baseload 
generation capacity available, the concept of storing energy becomes attractive, 
because energy storage can reduce the conventional peaking unit's (diesel and 
combustion turbines) requirements by allowing for the substitution of 
domestically available coal and uranium for the more expensive and Increasingly 
diff icult-to-obtain oil and natural gas. But storing electrical energy is both 
difficult and technologically limited. To date, pumped hydro has been the 
principal method of storing energy on a large scale; however, other methods 
of storing energy for utility application are being developed. One such method 
is thermal energy storage (TES) , utilizing the latent heat of phase-change 
storage media. 

The objectives of this task were to — 

• Identify a range of operation, 1 latent-heat TES temperatures required 
for application with fossil fuel-fired steam-electric plants 

• Define and identify a high-temperature TES concept that could utilize this 
high-temperature energy. 

These objectives were accomplished by evaluating peaking energy availability 
from coal- and nuclear-fueled baseload plants, design criteria and operating 
and economic requirements of the TES unit, and the general status of the 
technology of TES systems. In considering peaking TES subsystems, particular 
attention was given to assuring minimum Impact on the operating procedures of 
the primary generating plant, minimum degradation of the charging energy in 
the TES subsystem, and a realistically good overall system efficiency. 

System Concept Developmen t 

Several concepts exist for integrating TES peaking systems into power 
plants. For example,', a recent study by Grumman examined 10 different TES 
locations within nuclear and fossil fuel power plants, and Honeywell, Xerox, 
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and others* have proposed TES concepts for integration Kith solar-pov/ered 
plants. In selecting a TES utilization concept for anax/sls, current literature 
on TES systems was reviewed, and a meeting was held with representatives 
of a major utility In order to define the selection criteria and limits that 
would prohibit TES application.'^’^ The criteria and limits defined were 


• Significant Increases in the complexity of plant operations are 
undesirable. 

• The introduction of additional thermal/mechanical stress on the 
primary plant equipment is undesirable. 

• Retrofit of current feedwater TES concepts, and possibly all other 
TES concepts that assume retrofit, may not be feasible because of 
downtime requirements for system installation. 

• When the TES charging-steam circuit is such that at some point it 
passes through the plant's cooling tower, a portion of the energy 
rejected may have to be charged to the peaking system if fair 
energy-accounting practices are employed. 

• Latent-heat storage is well suited to the boiling operation. 

Based on these criteria, a TES peaking system concept that would allow 
considerable freedom in the specific design of the storage portion of the 
system, while taking advantage of the available information, was chosen. This 
TES peaking concept, shown in Figure 1, employs a separate power conversion 
loop (electrical generation subsystem) and a TES subsystem. During charging, 
steam coming directly from the primary plant's boiler is supplied to the TES 
subsystem, where it is condensed and returned as boiler feedwater. When 
discharging, the TES subsystem serves as a boiler, providing steam to the 
electrical generation subsystem ^ a conventional multistage condensing steam 
turbine. 


With the separate power conversion loop TES peaking system concept, 
additional thermal/mechanical stress is not introduced into the primary plant's 
generation section, and increased control complexity is limited to throttling 
the charging and discharging fluid. 

For this concept, the primary plant's steam generator determines, to 
a great extent, the TES subsystem’s operating conditions. The operating 
pressure determines the temperature at which the steam will condense and, hence, 
the temperature at which a major portion of the energy is available during 
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ciiarging. The output steam temperature Influences the discharging output 
temperature, and the feedwater return temperature establishes a target ’’ 
temperature for the condensate returned from the TES subsystem during charging. 
Typical primary plant steam generator characteristics are shown generically 
in Table 1 for the two fossil- and two nuclear-fueled plants. 

Table 1. TYPICAL POWER PLANT STEAM-GENERATOR CONDITIONS 

Conventional 


Output: Temperature, Steam, (°?) 518 (1000) 

Return Temperature, Water, (®F) 2^49 (^t80) 

Pressure, Pa (psla) 1.7 X 10'’ (2500) 

While it is desirable (and most likely possible) to combine a TES system 
with a pressurized water reactor (PWR) or a boiling water reactor (BWR) nuclear 
power plant, the low output steam temperature and pressure of just over 260‘’C 
(500"F) and 4.826 X 10^ Pa to 6.957 X 10® Pa (700 to 1009 psia) , respectively, 
prohibit the use of a separate power conversion loop TES peaking system. The 
conversion efficiency during discharge of the TES peaking system would be 
very low. 

In contrast, the fossil-fueled plants have boiler output temperature and 
ppressures of about 538°C (1000°F) and 1.655 X 10^ to 2.413 X 10^ Pa (2400 to 
3500 psia), permitting a considerably higher efficiency during discharge. 

Thus, we have assumed the TES utilization concept to be mated with a power plant 
that can provide superheated steam at 538°C (1000°F) and 1.655 X 10 Pa^(2400 psia) 
Additionally, for a TES system designed to accept superheated steam and discharge 
water in the charging mode, a staged or phased TES subsystem will be needed to 
minimize the reduction in temperatu) j availability from charging to discharging. 

Simplistlcally, the TES subsystem can be viewed as a subsystem with a 
charge/discharge mode. Both modes can then be broken down into three stages, 
with the discharge mode performing the reverse of the charge mode, as follows: 


Supercritical 

53B (1000) 

260 (510) 

1,3 X 10^ (3500) 


Nuclear 


PWK7 


266 (510) 

20 ^» (^* 00 ) 

2.5 X 10^' (690; 


mr 

288 (550) 

171 (340) 

3.6 X 10^' (965) 
A80020302 
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Stage 

1 

2 

3 


Charging Mode 
Desuperheating 
Condensing 
Condensate Cooling 


Discharging Mode 

Superheating 

Boiling 

Feedwater Heating 


Although the desuperheating/superheating and condensate cooling/feedwater 
heating stages could employ non-isothermal energy storage in a sensible storage 
material with temperature gradients in the direction of charging fluid flow, 
tl)e establisliment and maintenance of these gradients in storage and during 
discharge to ensure the system's proper operation is uncertain. Thus, we have 
assumed all TES subsystera stages to be nearly isothermal, employing phase 
change materials. 


In summary, the TES utilization concept selected consists of an electrical 
generation subsystem composed of a conventional multistage condensing steam 
turbine and a TES subsystem, with the subsystem employing at least three 
isothermal energy storage modules operating in series. Its only coiTnection 
to the primary plant, which provides steam at 538°C (1000°F) and 1.655 X 10^ Pa 
(2400 psia) , is a charging steam line and a feedwater return line. 


Operaring Conditions 

Base Cas e - 316°C (600°F) Out l et Steam 

In the TES subsystem, the charging steam is received from the primary 
plant's boiler at 538°C (1000°F) and 1.655 X 10^ Pa (2400 psia); at this pressure, 
the saturation temperature of the charging steam is about 350°C (662°F) , which 
defines an upper boundary for the operating temperature of the condensing/boiling 
stage of the TES subsystem. The lower boundary is simply the highest temperature 
and associated pressure that can be obtained during discharge, allowing for a 
reasonable driving temperature difference during charging and discharging. 

Assuming a temperature difference of about 50°C (90°F) between the PCM in the 
condensing/boiling stage and the charging and discharging fluid, the saturation 
temperature of the discharge steam was set at 247°C (477®F), resulting in a 
system pressure of 3.792 X 10^ Pa (550 psia) and a melt temperature for the 
PCM of about 297“C (567‘’F). Operating temperatures of the remaining stages 
were set somewhat arbitrarily after a discharge output temperature was selected, 
with reasonable temperatures for the feedwater return during charging and 
feedwater input during discharging being IcnoTO. 
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For a base-case, three-stage TES subsystem, the final temperature of the 
discharging fluid was conservatively set at 316°C (600°F), with feedwater 
input set at 178*C (352°F). Figure 2 conceptually shows the TES subsystem and 
its charging and discharging condition. With input steam conditions of 316®C (600“F) 
and 3.792 X 10® (550 psia) to the electrical generation subsystem and exhaust 
conditions of 6.895 X 10^ Pa (1 psia) and 5% moisture, a tentative heat balance was 
calculated for a plant producing 100 MW of net electric power output (Figure 3). 

A heat rate of 1.79 X 10^ J/kWlir (16,972 Btu/kWhr) was obtaiiied for a thermal- 
to-electric cycle efficiency of 20.1%. Table 2 provides additional cycle 
information. 

The total heat required in storage to achieve an 8-hour discharge period 
at 100 MW is 1.43 X 10^^ J (13.576 X 10^ Btu) . Discharging mass is 6.355 X 10® kg 
(14,010,552 lb) with a charging mass of 7.12 X 10^ kg (15,691,818 lb). 

Discharging flow rate is 7.94 X 10^ kg/hr (1,751,319 Ib/hr) for 8 hours. At 
a flow rate equal to the discharging rate, charging would require 8.96 hours. 

Additional Designs: 427°C (800°F) and 454°C (850°F) Outlet Steam 

T emperatures 

The base-case TES unit was assumed to provide steam at 316°C (600“F) and 
3.792 X 10*^ Pa (550 psia) via a multi-temperature storage unit. This TES unit 
was depicted as having nominally three isothermal sections to provide preheating, 
vaporization, and superheating to the discharging fluid. It is desirable, 
for efficient operation of the steam turbine, to increase the output temperature 
of the discharge fluid; this increase requires division into non-iso thermal 
superheating and preheating sections or the division of these two sections 
into severa] isothermal sections. 

Using the several isothermal sections approach, the number of isothermal 
storage sections and the charging-mass to discharging-mass ratio were varied 
to arrive at conceptual designs for TES units providing steam at 3.792 X 10® Pa 
(550 psia), 427°C (800°F) and 454°C (850°F) . Figures 4 and 5 show 
these designs. The 427 °C (800°F) case requires six isothermal sections 
and has a charge-to-discharge mass ratio of 1.2 while the 454"C (850®F) case 
requires seven sections with a charge-to-discharge mass ratio of 1.3. The 
increased charging-mass to discharging-mass ratios [the 316°C (600°F) unit’s 
1.12 ratio is the reference base] were required to obtain enough high-temperature 
energy in storage for superheating during discharge. 
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Table 2. HEAT-BALANCE INFORMATION FOR 100-MW PEAKING UNIT 
WITH 316 °C (600“F) OUTLET STEAM 


Full Throttle 

Turbine 


Steam Flow, kg/hr (Ib/hr) 


794,385 

(1,751,319) 

Pressure, Pa (psia) 


3.79 X 10® 

(550) 

Total Temperature, °C (®F) 


316 

(600) 

Condenser 

Steam Flow, kg/hr (Ib/hr) 


613,193 

(1,351,856) 

Pressure, Pa (psia) 


6894.8 

(1.0) 

3rd Extraction Heater 

Water Flow to Heater, kg/hr (Ib/hr) 


613,193 

(1,351,856) 

Turbi*'e Stage Pressure, Pa (psia) 


7.6 X 10^+ 

(11.0) 

Wati • Temperature Out, °C (°F) 


89 

(192) 

Water Temperature In, °C (°F) 
Extraction Enthalpy, J/kg (Btu/lb) 


39 

(102) 


2.7 X 10® 

(1160) 

Extraction Required, kg/hr (Ib/hr) 


55,695 

(122,786) 

2nd Extraction Heater 

Water Flow to Heater, kg/hr (Ib/hr) 


668,886 

(1,474,642) 

Turbine Stage Pressure, Pa (psia) 


3.17 X 10® 

(46) 

Water Temperature Out, “C (“F) 


133 

(272) 

Extraction Enthalpy, J/kg (Btu/lb) 


2.74 X 10® 

(1180) 

Extraction Required, kg/hr (Ib/hr) 


57,024 

(125.717) 

1st Extraction Heater 

Water Flow to Heater, kg/hr (Ib/hr) 


725,911 

(1,600,359) 

Turbine Stage Pressure, Pa (psia) 


1.0 X j.0® 

(145) 

Water Temperature Out, °C (“F) 


178' 

(352) 

Water Enthalpy In, J/kg (Btu/lb) 


7.54 X 10®. 

(324) 

Extraction Enthalpy, J/kg (Btu/lb) 


2.8 X 10® 

(1205) 

Extraction Required, kg/hr (Ib/hr) 


68,474 

(150,960) 

Boiler Feed, kg/hr (Ib/hr) 


79,438 

(1,751,319) 

Turbine Generator 

Enthalpy at Throttle, J/kg (Btu/lb) 


3.0 X 10® 

(1293) 

Work up to 1st Extraction Point, J/kg 

(Btu/lb) 

2.0 X 10® 

(88) 

Work up to 2nd Extraction Point, J/kg 

(Btu/lb) 

2.6 X 10® 

(113) 

Work up to 3rd Extraction Point, J/kg 

(Btu/lb) 

3.1 X 10® 

(133) 

Exhaust Enthalpy, J/kg (Btu/lb) 


2.4 X 10® 

(1053) 

Plant Performance 

Net Plant Sendout, kW 


100,000 

(100,000) 

Plant Heat Supplied, 10^ j/hr (10® Btu/hr) 

1.79 X 10® 

(1697.2) 

Overall Plant Heat Rate, J/kWhr (Btu/kWhr) 

1.79 X 10^ 

(16,972) 
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Tentative heat balances for a 100-MW steam turbine with three stages of 
water preheat were calculated for steam at 3.792 X 10® (550 psia) and at 
temperatures of 427" and 454"C (800" and 850"F) , The resulting system 
efficiencies were approximately 28% and 30%, a substantial Improvement over 
the 20% efficiency of the 316"C (600"F) base case. This Improvement is 
significant considering that a 30% efficient TES peaking system operating 
with a 38% efficient primary plant has an effective efficiency equal to that 
of pumped hydro with an electrlclty-to-electrlcity efficiency of about 70%. 

These heat balances and additional cycle Information are shown in Figure 6, 

Table 3, and Figure 7, Table 4, respectively. 

Our original thinking was that the TES subsystem designs having output 
temperatures of 316"C (600°F) and 454"C (850"F) would provide a conservative 
to somewhat optimistic parameterization. But this may not be the case because 
the lowest log-mean temperature difference for the 454"C (850"F) case is 
about 39"C (70°F) . Higher output temperatures could be obtained If lower 
log-mean temperature differences were affordable and higher charge-to-discharge 
mass ratios were acceptable. In addition, the system could be redesigned with 
more feedwater heating stages in the electrical generation subsystem, which 
would increase the efficiency of electrical generation and tend to eliminate 
part of the low-temperature sections of the TES subsystem. However, the 
temperature of the charging stream as it leaves the TES subsystem would 
Increase (possibly increasing the total plant operation complexity) to fully 
utilize the remaining energy in the stream. 

Economic Analysis 

Because cost data for 'he TES subsystem are uncertain, an equation 

defining output energy cost as a function of the peaking system's total annual 

cost and energy output was given to establish allowable costs for the TES 

subsystem. This analysis i-vis performed parametrically for a range of output 

energy costs ($30 to $80 milllon/kWhr ) , an electrical generation subsystem 

0 

cost of $200 and $400 per installed kilowatt, and levellzed fixed charge 
rates of 15% and 19% recommended by the Electric Power Research Institute (EPRI).^ 
These fixed charge rates represent tax preference cases for a plant life of 
20 years and no-tax preference cases, respectively. EPRI has recommended that 
the tax preference be excluded from long-term analyses, included for near-term 
analyses. Table 5 shows the EPRI-recommended levelized fixed charge rates 
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Table 1. IIKAT-BALANCIS INFORMATION FOR lOO-MW PEAKING UNIT 
WITH 427°C (800 “F) OUTLET STEAM 




Full 

Throttle 

Turbine 

Steam Flow, kg/hr (Ib/lir) 


50,847 

(1,120,989) 

PrcaHuro, Pa (psia) 


3.79 X 10*’ 

(550) 

Total Temperature, °C (°F) 


427 

(800) 

Condenser 

Steam Flow, kg/hr (ib/hr) 


39,246 

(865,227) 

Pressure, Pa (pvSia) 


6894.8 

(1.0) 

3rd Extraction Heater 

Water Flow to Heater, kg/hr (Ib/hr) 


39,246 

(865,227) 

Tiirbiue Stage Pressure, Pa (psia) 


7.6 X lO" 

(11.0) 

Water Temperature Out, °C (°F) 


89 

(192) 

Water Temperature In, °C (°F) 


39 

(102) 

Extraction Enthalpy, d/kg (Btu/lb) 


2.7 X 10 

(1160) 

Extraction Recjuirod , kg/hr (Ib/hr) 


35,682 

(78,665) 

2nd Ertractlon Heater 

Water Flow to Heater, kg/hr (Ib/hr) 


428,142 

(943,892) 

Turbine Stage Pressure, Pa (psia) 


3.17 X 10 

(46) 

Water Temperature Out, “C (°F) 


133 

(272) 

Extraction Enthalpy, J/kg (Btu/lb) 


2.74 X 10 

(1180) 

Extraction Required, kg/hr (Ib/hr) 


36,500 

(80,469) 

1st Extraction Heater 

Water Flow to Heater, kg/hr (Ib/hr) 


464,642 

(1,024,361) 

Turbine Stage Pressure, Pa (psia) 


1.0 X 10^ 

(145) 

Water Temperature Out, °C (°F) 


178 

7.54 X 10^ 

(352) 

Water Enthalpy In, d/kg (Btu/lb) 


(324) 

Extraction Enthalpy, J/kg (Btu/lb) 


2.8 X 10 

(1205) 

Extraction Required, kg/hr (Ib/hr) 


43,829 

(96,627) 

Boiler Feed, kg/hr (Ib/hr) 


508,472 

(1,120,989) 

Turbina Generator 

Enthalpy at Throttle, J/kg (Btu/lb) 


3. 28 X 10 

(1410) 

Work up to 1st Extraction Point, J/kg 

(Btu/lb) 

4.77 X 10 1 

(205) 

Work up to 2nd Extraction Point, J/kg 

(Btu/lb) 

5.00 X 10 

(215) 

Work up to 3rd Extraction Point, J/kg 

(Btu/lb) 

5.93 X 10 

(255) 

Exhaust Enthalpy, J/kg (Btu/lb) 


2.45 X lO*^’ 

(1053) 

Plant Performance 

Net Plant Sendout, kW 


100,000 

(100,000) 

Plant Heat Supplied, iO^'' j/hr (10^ Btu/hr) 

1.28 X 10*’ 

(1217.5) 

Overall Plant Heat Rate, j/kWhr (Btu/kWlir) 

1.28 X 10^ 

(12,175) 
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Table 4. HEAT- BALANCE INFORMATION FOR 100-MW PEAKING UNIT 
WITH 454“C (850“F) OUTLET STEAM 


Full Throttle 

Turbine 


Steam Flow, kg/hr (Ib/hr) 


465,926 

(1,027,192) 

Pressure, Pa (psla) 


3.79 X 10® 

(550) 

Total Temperature, °C ("F) 


454 

(850) 

Condenser 




Steam Flow, kg/hr (Ib/hr) 


359,652 

(792,897) 

Pressure, Pa (psia) 


6,894.8 

(1.0) 

3rd Extraction Heater 




Water Flow to Heater, kg/hr (Ib/hr) 


359,652 

(792,897) 

Turbine Stage Pressure, Pa (psia) 


7.6 X lO** 

(11.0) 

Water Temperature Out, ®C (®F) 


89 

(192.0) 

Water Temperature In, ®C (®F) 


39 

(102.0) 

Extraction Enthalpy, J/kg (Btu/lb) 


2.7 X 10^ 

(1160) 

Extraction Required, kg/hr (Ib/hr) 


32,666 

(72,017) 

2nd Extraction Heater 




Water Flow to Heater, kg/hr (Ib/hr) 


392,318 

(864,914) 

Turbine Stage Pressure, Pa (psia) 


3.17 X 10^ 

(46) 

Water Temperature Out, ®C (“F) 


133 

(272) 

Extraction Enthalpy, J/kg (Btu/lb) 


2.74 X 10^ 

(1180) 

Extraction Required, kg/hr (Ib/hr) 


33,446 

(73,736) 

1st Extraction Heater 




Water Flow to Heater, kg/hr (Ib/hr) 


425,764 

(938,650) 

Turbine Stage Pressure, Pa (psia) 


1.0 X 10® 

(145) 

Water Temperature Out, °C (“F) 


178 

(352) 

Water Enthalpy In, J/kg (Btu/lb) 


7.54 X 10® 

(324) 

Extraction Enthalpy, J/kg (Btu/lb 


2.8 X 10® 

(1205) 

Extraction Required, kg/hr (lb/’’r) 


40,162 

(88,542) 

Boiler Feed, kg/hr (Ib/hr) 


465,926 

(1,027,192) 

Turbine Generator 




Enthalpy at Throttle, J/kg (Btu/lb) 


3.35 X 10® 

(1440) 

Work up to 1st Extraction Point, J/kg 

(Btu/lb) 

5.40 X 10® 

(232) 

Work up to 2nd Extraction Point, J/kg 

(Btu/lb) 

5.70 X 10® 

(245) 

Work up to 3rd Extraction Point, J/kg 

(Btu/lb) 

6.63 X 10® 

(285) 

Exhaust Enthalpy, J/kg (Btu/lb) 


2.45 X 10® 

(1053) 

Plant Performance 




Net Plant Sendout, kW 


100,000 

(100,000) 

Plant Heat Supplied, 10^ J/hr (10^ Btu/hr) 

1.21 X 10® 

(1146.4) 

Overall Plant Heat Rate, J/kWhr (Btu/kWhr) 

1.21 X 10'^ 

(11,464) 
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a function of assumed plant life. Additional assumptions made In this analysis 
are as follows: 


• The peaking system size Is 100 MW^. 

• The peaking system Is designed for operation at full capacity 
for 8 hours per day. 

• The electrical generation subsystem investment cost Is Installed 
cost and Includes turbine/generator, condenser, and feedwater heaters. 

• The annual operating and maintenance cost for the electrical generator 
subsystem Is $0.003/kWhr. 

• The annual operating and maintenance cost for the TES subsystem Is 
3% of the capital Investment. 

• The TES subsystem energy losses are zero. 




Plant availability Is 85%. 


Table 5. 

EPRI -RECOMMENDED LEVELIZED FIXED CHARGE RATES ^ 

Facility Life, 

Levelized 

Fixed Charges , %* 

Years 

No Tax Preference 

With Tax Preference 

5 

34 

29 

10 

23 

19 

15 

20 

16 

20 

19 

15 

25 

18 

15 

30 

18 

15 

35 

18 

15 

40 

18 

15 

45 

18 

15 

50 

18 

15 


Includes return, depreciation, allowance for retirement dispersion, income 
taxes, other taxes, and insurance. Excludes operation and maintenance. 

Essentially this economic analysis treated the break-even cost of such 
a system as an unknown and solved for it in terms of the electrical generation 
subsystem and desired electric power costs. 
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The derivation is as follows: 



EC 

!= 

AC 

S X. PF 


(1) 


AC 

- 

ACT + ACE 


(2) 

ACE = 

(FD X 

IE) 

+ (FN X WC) 

+ AE 

(3) 


IE = 

El 

X S X io“^ 


(4) 

WC = 

0.02 X 

IE 

(See assumptions.) 

(5) 

(3.0 

X S X 

PF 

X 8760) lo“^ 

(See assumptions.) 

(6) 


where 

EC energy cost, $ inlllion/kWhr 

K = conversion constant 

S = plant size, MIV 

PF ~ plant factor (decimal) 

AC = annual plant cost, $ million 

ACT “ annual TES subsystem cost, $ million 

ACE = annual electrical generation subsystem cost, $ million 
FD = fixed charge rate for depreciating investment 

FN = fixed charge rate for nondepreciating investment 

WC = working capital 

AE =5 annual operating and maintenance cost for electrical 
generation subsystem, $ million 

IE = total investment for electrical generation subsystem, 

$ million 

El » electrical generation subsystem cost, $/kW 

combining 

ACT = g . A EC _ Jg) ^ ^ ^ 

and substituting values 

ACT = S jPF[ (8.756 X 10"3 X EC) - 2.628 X lO"^] _[i . 8235 X lO"'* X El ]| (8) 
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Assuming the same fixed charge rates and percentage of Investment for 
working capital as before, with maintenance on the TES subsystem taken at 3 % 
of the capital Investment, the break-eVen cost per kilowatt : electric for the 
TES subsystem can be stated as — 


ET = 


ACT - FC 
(2.1235 X 10"^)S 


(9) 


or 


» act - FC 
2.1235 X 10”^ 


( 10 ) 


where 

ET * break-even thermal energy storage cost, $/kWg 
ACT = break-even annual cost of the subsystem, $ million 
S = system size, MW 

FC = annual fuel cost for charging TES unit, $ million 


The results are shown in Figures 8, 9, and 10 for the 316“C (600®F) , 

427“C (800°F) and 454°C (850°F) systems. The left axis shows the total annual 
break-even cost of the TES subsystem including charging fuel costs. The 
right axis shows the total break-even capital investment for the TES subsystem 
per installed kilowatt-hour (electric) for a charging fuel cost of j 

($1/10^ Btu) (total charging energy costs). The zero point of the right axis 
is the point where the allowed annual TES subsystem cost equals the annual 
fuel cost. The investment costs shown on this axis are equal to $6.94 and 
$5.68 per million dollars of allowed annual cost above the annual fuel cost 
for the 15% and 19% fixed charged rates, respectively. 

For the 454®C (850“F) case, the break-even storage cost with electrical 
generation subsystem cost of $200/kW and desired energy costs of $60 million/kWhrg 
are $42/kV7hre and $57/kWhrg (19% and 15% fixed charged rate). With the 316®C 
(600“F) case, these costs are reduced to $36/kWhrg and $48/kWhrg. This 
$8/kWhrg to $9/kWhre reduction in break-even storage cost for the 316®C (600“F) 
case results from the $1.4 X 10® increase in annual fuel costs, a result 
of lower efficiency. 
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Figure 8. TARGET COSTS FOR A TES SUBSYSTEM AS A FUNCTION OF ELECTRICAL 
GENERATION SUBSYSTEM COST, FIXED CHARGE RATE, AND OUTPUT ELECTRICITY 
COST WITH "^le^C (600'’F) OUTLET STEAM 
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Figure 9. TARGET COSTS FOR A TES SUBSYSTEM AS A FUNCTION OF ELECTRICAL 
GENERATION SUBSYSTEM COST, FIXED CHARGE RATE, AND OUTPUT ELECTRICITY 
COST WITH 427 “C (800°F) OUTLET STEAM 
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Status of Technology 

Numerous programs to study the applicability of using the latent heat 
properties of various inorganic salts to store thermal energy in the respective 
phase-change temperature regimes are either complete or underway. 

Storage of thermal energy in PCM is desirable because large quantities of 
energy are involved in the phase-change process, and this energy is transferred 
(ideally) at a consf'.ant temperature (the phase-change or melting point). As 
a result of these programs, salts and salt mixtures have been identified for 
applications spanning a temperature range of 93‘’C (200°F) to in excess of 
1090°C (2000°F). The lower range of this temperature span includes space 
heating applications, and the upper temperature range includes several solar 
and vehicular propulsion applications. Storage of thermal energy from cyclic 
industrial processes and in electric utility steam power plants falls between 
the two temperature extremes. Each application in this wide temperature range 
has its own unique set of requirements. For example, in the application where 
thermal energy is stored and subsequently used to power a vehicle, weight 
and volume of the storage subsystem is of prime importance. However, in* 
identifying salts for utility applications, the prime factors are cost, 
reliability, safety, and useful operating lifetimes. 

The major portion of the TES subsystem consists of the TES salt and the 
containment vessel/HX. These two factors must be considered simultaneously 
because of their various interactions. For example, a salt chosen for 
its low cost and high heat of fusion may be corrosive toward low 
cost containment vessel/HX materials, necessitating thO use of more expensive 
and/or complex materials (e.g., linings) or the selection of a different salt. 
Similarly, if the chosen salt readily reacts with the environment (air or H 2 O) , 
then special containment and sealing arrangements must be provided, 

Experimental work at several laboratories (IGT, Comstock and Westcott, 
Minneapolis Honeywell, and Philips) ^ has addresBed these trade-offs; 
however, insufficient data are available to provide exact and quantitative 
numbers, especially long-term salt/containment/HX data. 

The design and costs of the HX required in TES systems utilizing PCM are 
governed by the complexity of the heat transfer processes occurring as well as 
the corrosivity of the salt. The latter factor can be rather easily accounted 
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for by materials compatibility testing. But the former presents fundamental 
and technical constraints that could greatly Increase the HX complexity. This 
complexity results from the fact that as thermal energy is removed from the 
molten salt, a layer of solidified salt becomes deposited on the HX surface. 
Because of the poor thermal conducltivity of the solidified salt, the heat 
exchange rate (heat flux) decreases with time as the layer thickness increases. 
Therefore, although ideally the heat flux from an isothermal source should 
remain constant, the observed heat flux from a phase-change TES system varies, 
approaching that more typical of a sensible heat storage source. Programs 
are now under way to model, predict, and experimentally verify the rates of 
heat transfer in TES PCM.^°>’^ The models developed so far for selected 
conditions predict quite accurately the observed experimental behavior as the 
salt layer thickness increases. Currently, most experimental efforts focus 
on ways to minimize the salt layer thickness or improve its thermal 
conductivity. Physical as well as chemical means to control the salt layer 
thickness are being utilized. Some physical means utilize scrapers and others 
use vibratory techniques to remove the solidified layer. Modifying the 
chemical composition of the salt, so that a "slush" and weakly crystallized 
layer rather than a solidified layer is formed, is also being tested. Because 
the thermal conductivity of most salts is only of the order of 1.73 W/m-K 
(1 Btu/hr-ft-“F) , the addition of a more conductive material having a cellular 
or directional conductivity is being investigated. Aluminum, for example, is 
a preferred conductivity promoter because it is approximately two orders of 
magnitude higher in thermal conductivity. Because of the state of development 
of some of these advanced heat exchange concepts, comparative quantitative heat 
flux data are not available. 

The technological readiness of passive HX concepts for the TES subsystem 
discussed in this paper can be assessed if one makes the following simplifying 
assumption and calculates parametrically the cost breakdown for the TES sub- 
system components — the containment vessel/HX and TES salt. 

• Operational concerns, such as reliability, safety, and useful life, are 

solvable; these costs are included in the containment vessel/HX. cost. 
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Tn performing this calculation, we selected the 454®C (850®F) case as the 
reference design and $60 million/kWhr0 as the target energy cost. Thus, the 
break-even installed cost of the TES subsystem is $52/kWhr (average of the 
15% and 19% fixed-charge rate break-even costs, with electrical generation sub- 
system costing $200/kW). 

Further assumptions were — 

• Installed cost is equal to 150% of materials cost. Thus break-even 
materials cost for 1 kWhr of electric storage is about $34.60. 

• Passive heat exchange designs have heat flux rates of 3152 W/m^ 

(1000 Btu/hr-ft2). 

• The TES salt has a heat fusion of 2.326 X 10^ j/kg (100 Btu/lb) , 
with a cost of $0. 44/kg ($0. 20/lb). 

The result of this analysis shows that for passive heat exchangers having 
a heat fl.ix of 3152 W/m^ (1000 Btu/hr-ft^), the cost ceiling for containment 
vessel/HX is about $86/m^ ($8/ft^) of heat exchange area. Expressing these 
break-even costs in terms of heat flux, the ceiling cost is $0.26/100 W/m^ 
($0,817/100 Btu/hr-ft^). Further analysis using different salt costs showed 
that a $0. 11/kg ($0. 05/lb) change in salt cost results in $0.13/100 W/m^ 
($0.40/100 Btu/hr-ft^) change in the ceiling cost of the containment vessel/HX. 
Table 6 shows containment vessel/HX ceiling costs for several salt costs. 
Overall subsystem costs may be either salt-cost intensive or HX-cost intensive, 
depending on the application's duty requirements, which dictate selection of 
salt properties and HX configuration. 
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Table 6. CONTAINMENT VESSEL/HX CEILING COSTS AS A FUNCTION OF SALT COST* 
(Expressed in Terms of HX Heat Flux) 


Salt Cost Containment Vessel/HX Cost 


$/kg 

($/lb)** 

J/d 

(Btu/C)*" 

$/100 W/m2 

($/100 Btu/hr-ft^) 

0.22 

(0.10) 

10,551 

(10) 

0.513 

(1.617) 

0.33 

(0.15) 

7,037 

(6.67) 

0.386 

(1.217) 

0.44 

(0.20) 

5,275 

(5) 

0.249 

(0.817) 

0.55 

(0.25) 

4,220 

(4) 

0.132 

(0.017) 


Based on a thermal-to-electric conversion efficiency of 29.8% and an 
allowed storage cost of $52/kWhrg installed ($34.6/kWhrg not Installed). 

Assumes salt heat of fusion of 2.326 X 10^ J/kg (100 Btu/lb). 

Expressing a salt's cost and heat of fusion in this manner permits any 
combination of salt cost and !.aat of fusion to be used in this table 
[e.g., $0. 22/kg ($0.1/lb) with 2.326 X 10^ J/kg (100 Btu/lb) - 10,573 
J/(? (10 Btu/<?); $0. 33/kg ($0. 15/lb) with 3.489 X 10^ J/kg (150 Btu/lb) = 
10,573 J/O (10 Btu/<?)]. 

Summary 

TES peaking systems for utility applications employing a separate power 
conversion loop appear viable when interfaced with high-temperature, high- 
pressure base-load steam-electric plants. Their Impact on the primary plant's 
operation is minimal, but storage, efficiency is high. A 30% efficient TES 
peaking system operating with a 38% efficient primary plant has an effective 
efficiency equal to that of pumped hydro with an electricity-to-electricity 
efficiency of about 70%. Nevertheless, with current passive HX costs 
in the $215 to $269/m2 ($20 to $25/ft2) range, the economic viability of the 
TES peaking system clearly rests in the development of reasonably priced 
containment vessels/heat exchangers with high heat flux rates. 
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Task 2, TES Modi;le Testing 

2.1 Engineering-Scale Testing 
Introduction 

Under IGT's previous thermal energy storage development contract 
(ERDA EY-76-C-02-2888, "Molten Salt Thermal Energy Storage Systems") three 
8 kybrj-h'^^apacity engineering -scale TES modules were successfully de- 
signed, fabricated, a... operated with LIKCO3 PCM with a mid-range 
melting point of 505®C (941°F),^‘* One of these units (Unit 2) was operated 
for 1400 hours and 30 thermal cycles as part of that effort. The operation 
of Unit 2 wns voluntarily suspended on June 3, 1977 and it was partially 
disassembled to evaluate the condition of the system. The containment, HX, 
and salt were found to be in excellent condition. Downtime was approximately 
five months., during which the solidified LiKC0;3 salt was in contact with 
ambient air with no detrimental effects. Based on the promising thermal 
performance results and short-term stability of this TES system, the unit 
was reassembled and operated to obtain long-term engineering and compatibility 
data under the present program. 

The specific objectives of this subtask were to — 

• Establish the endurance and performance of this engineering-scale 
TES system during prolonged charge/discharge cycling 

• Determine the effc:_ts of long-term, intermediate temperature 
exposure and cyclic temperature operation on the phase-change, 
salt, containment, and HX materials. 

The objectives were met after successfully operating Unit 2 for 5650 hours and 
129 thermal cycles at 480° to 535° (896° to 995°F) . Very good chemical compati- 
bility was demonstrated between the LiKC03 salt and the austenitic Type 316 
stainless-steel (SS) materials of construction, with the carbonate salt exhibiting 
no hysteresis effects. 

Cyclic Endurance Testing 

The engineering-scale TES containment design, test-station module, and 
overall system schematic are shown in Figures 11, 12, and 13. The unit 
utilizes a single-pass, single-tube HX design and contains 68 kg (150 lb) of 
reagent-grade LiKC03 salt material (35 wt % L12C03 - 65 wt % K2CO3). This 
salt has a melting point of 505°C (941°F) and a heat of fusion of 3.44 X 10^ J/kg 
(148 Btu/lb). More detailed design and fabrication data were presented in a 
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Figure 11. ENGINEERING- SCALE THERMAL ENERGY STORAGE UNIT 
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previous report. The system was thermally charged by a resistance heater 
located In the HX tube until the salt was molten and uniformly heated to 
535“C (995®F) , The system was discharged when a hlgh-veloclty stream of ambient 
air at a mass flow rate of almost 159 kg/hr (350 Ib/hr) was passed up 
through the HX tube. The LIKCO3 PCM, which was at a temperature of 
535®C (995“F) Initially, was thus cooled through the 505®C (941“F) solidification 
point to a temperature of approximately 480“C (896“F) . The unit was then 
thermally recharged overnight In preparation for another run. 

Following the 5 1/2-month shutdown from June to November 1977, the 
endurance-testing phase for Unit No. 2 was continued on November 16, 1977. 

Cyclic endurance testing was continued during the following 6 months, when 
Unit No. 2 accumulated a total of 129 charge/discharge cycles and 5650 hours 
of operation. The unit was terminated on May 12, 1978. 


The discharge performance of Unit 2 was evaluated quantitatively 
periodically to establish the long-term cycling stability of a molten carbonate 
TES module. The heat flux, Q, from the solld/llquld Interface In the salt to 
the air working fluid, was determined as a function of discharge time by 
monitoring the Inlet and outlet air temperatures and calculating Q from the 
expression — 

mCpAT 

Q - ~~i — (11) 

where Q = heat flux, W/m^ (Btu/hr-ft^) 

m = air mass flow rate, kg/hr (Ib/hr) 


Cp = heat capacity of air, J/kg-®C (Btu/lb-®F) 

AT = air-temperature differential = T^^^ — T^^, ®C (°F) 
A = heat-transfer area, m^ (ft^)^** 


Results of such analyses are presented In Figure 14, which shows the 
variation of heat flux with discharge time for five cycles at different points 
during the lifetime of this unit. Table 7 summarizes the time required to 
discharge the system and the time-averaged heat flux determined for each of 
these cycles. Discharge times ranged from 315 to 345 minutes, with time-averaged 
heat fluxes ranging from 1,17 to 1.31 X 10** W/m^ (3715 to 4160 Btu/hr-ft^). 
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These data define a rather narrow discharge performance band, indicating no 
significant loss of thermal energy-storage capacity or heat-transfer capability 
during more than 5000 hours of operation. This is a significant observation 
because it indicates dependable performance does occur after a prolonged 
shutdown, which would occur in a utility application for equlment maintenance. 

Table 7. EXPERIMENTAL RESULTS FOR DISCHARGING OF ENGINEERING-SCALE TES UNIT 2 

(LIKCO3 Salt) 


Discharge 

Cycle 

Number 

Cumulative 
Operating 
Time, hr 

Time Required 
for Discharge, 
min 

Average Heat Flux 
W/m^ (Btu/hr~ft^) 

2-15 

500 

315 

13,117 

(4161) 

2-65 

2900 

330 

12,373 

(3925) 

2-81 

3550 

330 

10,570 

(3353) 

2-102 

4400 

330 

11,711 

(3715) 

2-127 

5360 

345 

11,831 

(3753) 


Post-Test Examination 

After termination. Unit 2 was dismantled to evaluate the stability of 
the LIKCO3 salt and the nature and extent of corrosion of the stainless-steel 
HX tube and containment vessel. A side-view photograph of the 0.3 m (12 in.) 
in diameter X 0.46 m (18 in.) -high Unit 2 module after termination and removal 
of the insulation is shown in Figure 15. The appearance of the upper surface 
of solidified LIKCO3 salt after removal of the top cover is shown in Figure 16. 
The carbonate salt was very clean, free of any significant discoloration from 
chemical interactions with the stainless-steel HX tube and containment vessel. 

A very thin and weakly adherent film of corrosion product was observed 
on the inner wall of the containment vessel at tlie salt/gas interface. The 
film had become detached from some areas and formed a few small localized 
reddish-brown deposits on the surface of the solid salt. Small samples of 
the corrosion-product scale were scraped from the containment surface and 
washed with an acetic acld/acetlc anhydride solution to remove the soluble 
salts. The residue was collected on a millipore filter paper and analyzed by 
X-ray dlf fractometry (Cu-K^ radiation) in an attempt to identify the compound(s) 
present. As shown in Table 8, six distinct peaks were observed having d-spacings 
ranging from 4.23 to 2.05 angstroms. Comparison of these peaks with 
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Figure 16. APPEARANCE OF LiKCOa SALT IN ENGINEERING 
UNIT 2 AFTER ] 29 CYCLES AND 5650 HOURS OF ELEVATED- 
TEMPERATURE OPERATION 
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available literature data for lithium-containing mixed oxide compounds indicates 
that LlCr02 or LiFe02 may have formed on the Type 316 SS. This is a 
plausible finding based on the nominal composition of this steel (17% Cr, 

12% Nl, 2.5% Mo, 2% Mn, 1% Si, balance Fe, by weight) although the amount of 
sample available was too small to obtain a high-precision X-ray diffraction 
pattern or confirmatory chemical analyses. 


Table 8. X-W DIFFRACTION ANALYSIS OF CORROSION-PRODUCT SCALE FROM 

ENGINEERING UNIT 2 


Corrosion- 

-Product 

a-LiFeO, (17-938)’' 

8-LiFeO, 

(17-936)* 

LiCrO^ 

(24-601) 

d(A) 

I/I, 

d(A) 

I/Io 

d(A) 

l/lo 

d(A) 

I/Io 

4.23 

45 





4.800 

75 

3.81 

13 







^.72 

8 







2.42 

4 

2.394 

50 

2.392 

30 

2.473 

45 

2.09 

13 



2.136 

70 



2.05 

100 

2.073 

100 

2.043 

100 

2.060 

100 


* Diffraction data for LlFe02 and LiCr02 compounds taken from Reference 15. 


The aluminum foil gasket used to form a non-hermetic seal between the 
interior of the vessel and the ambient atmosphere was in reasonably good 
condition, having experienced some oxidation, embrittlement, and cracking. 
(See Figure 16.) Visual examination of the module exterior after removal of 
the insulation revealed a small deposit of carbonate that was apparently 
transporated by the vapor phase from the melt through this aluminum gasket 
seal area between the cover and vessel. Visual examination also revealed 
that a small amount of carbonate was transported (most likely via the vapor 
phase) from the melt into the upper portion of the HX pipe, resulting in some 
minor corrosion on the inner surface of this pipe. A very slight leak was 
observed at one of the thermocouple fittings on the bottom of the vessel. 

But, in general, visual inspection indicated that the molten carbonate TES 
module was in very good condition after 5650 hours of elevated-temperature 
service. 
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Metallographlc analysis of the 316 stainless-steel vessel revealed that 
a stable oxide scale approximately 25pm (0.001 in.) thick was formed on the 
inner surface during the 5650 hours of exposure to the carbonate environment. 

An optical micrograph of the containment inner surface exposed to LiKCOs at 
mid-salt level is shown in Figure 17. In addition to the outer oxide layer, 
a sub-scale zone of Intergranular penetration approximately 76pm (0.003 in.) 
in e..cent was also observed. There was no visible evidence of spalling of the 
outer oxide scale as a result of repeated thermal cycling; the scale appears to 
have been providing protection against further oxidation of the stainless 
steel. No significant thermal aging effects were observed in the bulk stainless- 
steel material. 

Samples of the L 1 KC 03 salt were taken from various locations in the TES 
unit for analysis. Differential scanning calorimetry (DSC) was utilized to 
determine whether any changes in the thermal behavior of the LIKCO 3 salt; 
had occurred as a result of the prolonged operation and cycling at elevated 
temperature. The melting point of the salt after 5650 hours of testing was 
found by DSC to consistently be 502. 5®C (936. 5°F) compared to 505®C (941“F) 
for the unexposed material, indicating very good salt stability. The chemical 
analysis of these samples, appearing in Table 9, revealed that there was no 
change in the composition of the salt due to partial decomposition or chemical 
reaction with the HX or containment materials. 

Table 9. CHEMICAL, ANALYSIS OF LiKCOs SALT FROM ENGINEERING-SCALE UNIT 2 


COMPOSITION. Wt % 



Li 

K 

CO 3 

Fe 

Ni 

Cr 

PRE-TEST 

6.57 

36.78 

56.64 

<0.01 

<0.01 

<0.01 

POST-TEST 







2.54cm (1 in.)* 

6.58 

38.76 

55.27 

0.01 

<0.02 

<0.01 

5.08cm (2 in.)* 

6.43 

37.43 

57.01 

0.01 

<0.01 

<0.01 

6.35cm (2.5 in.)* 

6.72 

37.54 

57.48 

0.01 

<0.01 

<0.01 

Average Post-Test 

6.58 

37.91 

56.58 

0.01 

<0.01 

<0.01 

* 

Samples taken at mid-salt 

level , 2 . 

54, 5.08, 

and 6.35 

cm 




(1, 2, and 2.5 in.) away from containment OD. 
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Figure 17. TYPE 316 SS CONTAINMENT FROM ENGINEERING UNIT 2 AFTER 5650 HOURS 
IF TESTING, SALT SIDE, MID-HEIGHT (FeCl2-HCl-H20 Etchant) 
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The high degree of stability of this alkali carbonate TES module after 
5650 hours of operation is very encouraging > especially since this performance 
was achieved without a hermetic seal betv?een the ambient and gas phase above 
the PCM salt and after a 5-month exposure of the salt to ambient during 
a mid-test shutdown. 

2.2 Laboratory-Scale Testing 
Introduction 

The results of the system studies performed under Task 1 Indicated that the 
PCM for application in a steam-Ranklne cycle should have a transformation 
temperature lass than 350“C (662®F) . This temperature is defined by the 
saturation temperature of the steam working fluid at a pressure of 1.6 X 10^ Pa 
(2,400 psia) In the possible TES schemes indentlfied in Task 1, for example, 

67 to 77% of the total heat transferred between the PCM and the working fluid 
would be associated with a unit operating at an approximate temperature of 
SOO^C (570®F) . Because this temperature is lower than the lowest melting 
carbonate eutectic known [43.5 mol % Li 2 C 03 / 31.5 mol % Na 2 C 03 / 25.0 mol % 

K 2 CO 3 ternary; 397°C melting point (mp) ] , identification of a suitable PCM 
for this application would require the investigation of ways to suppress 
the ternary's mp, possibly by doping it with other carbonates (MgC 03 , 

SrC 03 , CaC 03 , BaC 03 ), Consideration of non-carbonate salts in this temperature 
regime was outside the scope of the present progi'.?;n. 

In view of the aforementioned considerations, the major experimental effort 
was redirected from the TES application originally intended (steam-Rankine 
peaking cycle) to solar-thermal systems applications in the temperature range 
538° to 870°C (1000° to 1600°F) . Carbonate salts are particuldrly attractive 
for such high temperature applications because of their relatively high heats 
of fusion, thermal stability and compatibility with materials of construction. 
Emphasis was given to salt compositions based on Li 2 C 03 , Na 2 C 03 , and K 2 CO 3 , 
from which congruently melting carbonate mixtures are available covering the 
temperature range 397°C (747°F) for the (,Li-Na-K) 2 C 03 ternary to 898°C (1648°F) 
for pure K 2 CO 3 . 
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Salt Selection 

The most important salt properties to be considered when selecting a 
thermal energy storage medium were Identified and discussed in detail 
previcasly. These factors Include — 

• Melting temperature 

• Heat of fusion 

• Heat capacity (solid and liquid) 

• Tliei'inal conductivity (solid and liquid) 

• Density 

• Volume change on fusion 

• Volumetric thermal expansion 

• Cost and availability 

• Chemical stability 

• Containment requirements 

• Safety 

The most reliable thermophysical and thermochemical carbonate property data 
available are for the pure alkali carbonates, L12C03, Na2C03, and K2CO3. 

These data have been reviewed by Janz,^^ Stern and Weise,^*^ and most recently 
by Janz et al.^^ Table 10 provides a summary of pertinent data for these 
pure carbonates. 

The heats of fusion of L12C03 (mp 723°C, 1333*^), Na^.COs (mp SSS^C, 
1576“F), and K2CO3 (mp 898®C, 1648'’P) were measured by Janz, Neuenschwander , 
and Kelly20 to be 6.07 X 10^ 2.65 X 10^ and 2.00 X lO^^ J/kg (261, 114, 
and 86 Btu/lb), respectively, using a drop calorimetric technique. Based 
upon these results and complementary calorimetry work by Rolin and Recapet, 
Janz et al.^^ most recently recommended the following heat of fusion values: 
LI2CO3 - 5.56 X 10^ J/kg (239 Btu/lb); Na2C03 - 2.81 X lO^ j/kg (121 Btu/lb); 
K2CO3 - 2.00 X 10® J/kg (86 Btu/lb). However, the former values were used 
during the course of this work. 
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Table 10. SELECTED PROPERTIES OF ALKALI CARBONATES 



LI 2 CO 3 

Na 2 C 03 

K 2 CO 3 

Molecular Weight 

73.89 

105.99 

138.21 

Melting Point, °C 

723 

858 

898 

Heat of Fusion, J/kg (Btu/lb) 

6.07 X 10^ 
(261) 

2.65 X 10^ 
(114) 

2.00 X 10^ 
( 86 ) 

Heat Capacity, J/mol-K 
(50®C above mp) 

188.1 

194.8 

208.6 

Density, kg/m^ (Ib/ft^) 
Solid at 25°C 
I Lquid at 50“C above rop 

2114 (132) 
1810 (113) 

2531 (158) 
1954 (122) 

2434 (152) 
1874 (117) 

Volume Change on Melting, % 

+6.9 

+16.2 

+16.4 

Thermal Conductivity, W/m-K 
(Btu/hr-ft-°P) 

Solid at 50*C below mp 
Liquid at 50°C above mp 

1.45 (2.51) 
2.14 (3.70) 

1.27 (2.19) 
1.88 (3.26) 

— 

Dissociation Constant, Pa (atm) 
(SO^C above mp) 

^ 2.6 

(4.2 X 10 ) 

0.06 „ 
(5.9 X 10 ” ) 

0.006 

(6.4 X lO"' 

Viscosity, Pa-s (cP) 
(50*C above mp) 

4.65 X 10 "^ 
(4.65) 

2.8 X 10 “^ 
(2.80) 

2.22 X 10 “^ 
( 2 . 22 ) 

Surface Tension, N/m 
(50°C above mp) 

0.242 

0.209 

0.166 

Technical-Grade Salt Coat. 
$/kg ($/lb) 

2.05 

(0.93) 

0.07 

(0.03) 

0.42 

(0.19) 
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Expressions for the heat capacities of the molten salts have been 

presented by Janz et al,^^ with Cp units of J/mol-K. 

_2 

LiaCOa Cp = 129.0 + 5.66 X 10 T(K) (12) 

NaaCOg Cp = 142.1 + 4.47 X 10“^‘ T(K) (13) 

KaCOa Cp = 154.5 + 4.45 X lO"^ T(K) (14) 

The enthalpy curves shown in Figure 18 represent the sensible and latent heat 
contents for the pure carbonates and estimates for several carbonate mixtures 
of interest, based on molar additivity. The sensible heat content curve of 
MgO is also shown for comparison. 

Thermal conductivities of solid and liquid LiaCOa and NaaCOa have been 
measured by Egorov and Revyakina^^ using the absolute steady-state radial 
heat flux method. Conductivity equations have been presented by Janz et al,^^ 
with units of k in W/m-K. 

LiaCOa k(s) = 7.59 -1.29 X 10“^ T(K) + 6.81 X 10~^ T(K)^ (700-1000K) (15) 

k(l) “ 1.36 X 10"3 + 3.33 X 10~3 T(K) (16) 

NaaCOa k(s) = 8.53 -1.40 X 10“^ T(K) + 6.75 X lO"^ T(K)^ (700-1135K) (17) 
k(l) = 4.94 X 10"^ + 1.18 X 10"^ T(K) (18) 

No thermal conductivity data have been reported for KaCOa. 

The densities of carbonate salts at 25®C (77°F) are 2110 kg/m^ for 
LiaCOa, 2530 kg/m^ for NaaCOa, and 2430 kg/m^ for KaCOa. Densities of 
the molten salts at a temperature 50°C (90°F) above the melting point are: 

1810 kg/m3 for LiaCOa at 773°C (1423°F) , 1950 kg/m^ for NaaCOa at 908°C 
(1666®F) , and 1870 kg/m^ for KaCOa at 948°C (1738“F). The volume changes on 
melting have been estimated to be 6.9%, 16.2%, and 16.4% for LiaCOa, NaaCOa, 
and KaCOa, respectively. 


An English unit to SI unit conversion table is presented in Appendix A. 
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The chemical stability of molten carbonates depends primarily on 
temperature and the gaseous environment maintained over the melt. Decomposition 
of the molten metal carbonate M 2 C 03 (Jl) can be expressed as ~ 

M2C03(li) t M20(«) + C02(g). (19) 

The thermodynamic equilibrium between the melt and the gas phase over it is 

described In terms of the dissociation constant, K. ~ 

d 

^C02 (e) 

K = (20) 

^I 2 C 03 (O 

The term in the above equation is the equilibrium GO 2 dissociation pressure, 

and ^j 20 (<) ^InC 03 ('!) thermodynamic activities of tl>e metal oxide 

and metal carbonate, respectively. For ideal solutions, the activity terms 
may be replaced by the mole fraction (X) of each species. Dissociation constants 
for L 1 ;iC 03 , Na 2 C 03 , and K 2 CO 3 at 50“C (90°F) above their melting points have 
been estimated as 42.6 Pa, 0.06 Fa, and 0.006 Pa (4.2 X 10 ** atm, 5.9 X 10~^ 
atm, and 6.4 X 10 ^ atm), respectively. Thus, for a given oxide content and 
constant degree of superheat above the melting point, the stability is in the 
order — 

K2CO3 > Na^COs > Li;<C03. 

Partial decomposition of molten carbonates with the formation of the oxide 

ion species, which is more chemically ayyreasive toward staii^less steel than 

carbonate ion, can be minimized by maintenance of a finite GO?, partial pressure 

-3 

over the melt. An M 2 O mole fraction of 10 would be established in molten 
hi. 2 C 03 , Na 2 C 03 , and K 2 GO 3 under GO 2 pai'tial pressures of 4.3 X 10*^ Pa, 59.8 Pa, 
and 6.48 Pa (0.42 atm, 5.9 X 10 atm, and 6.4 X 10 atm), respectively, at 
50°C (90°F) above the melting points. No systematic studies of vaporization 
have been reported for Li 2 G 03 , Na 2 C 0 a or K 2 GO 3 as carbonate - 

MjCOgCO ^ M 2 G 03 (g) (21 ; 
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Expressions have also been presented by Janz et al^^ for the viscosity 
(n, Pa-s) and surface tension (y» N/m) of the molten alkali carbonates -• 

LiaCOa: n » 5.259 + 0.015 T(K) - 13.858 X 10~^ T(K)2 + 4.313 X lO"^ T(K)^(22) 


y « 0.285 - 4.06 X lO"^ T(K) (23) 

Na-COa: n “ 3.832 X lO"^ exp [26,260/RT(K) ] (24) 

Y = 0.269 - 5.02 X lO"^ T(K) (25) 

K 2 CO 3 : n = 1.161 X lO"® exp t29,487/RT(K)] (26) 

Y = 0.338 - 21.71 X lo"^ T(K) + 6.25 X 10‘® T(K)^ (27) 


Selected values of these properties are summarized in Table 10. 

The alkaline earth carbonates MgCOa, CaCOa , SrCOa, and BaCOa are of 
interest as additives to alkali carbonate PCM mixtures because of their 
ability to modify the melting and solidification behavior. The pure alkaline 
earth carbonates are generally less stable than the alkali carbonates against 
thermal decomposition, as indicated by their dissociation constants given in 
Table 11. As pure compounds they generally decompose to M 2 O and CO 2 rather 
than melting. Nevertheless, their decomposition in molten mixtures with 
alkali carbonates can be adequately controlled in certain cases by providing 
a CO 2 partial pressure over the melts. 

Table 11. SELECTED PROPERTIES OF ALKALINE EARTH CARBONATES 





I 





MgCOa 

cacoa 

SrCOa 

BaCOa 

Molecular Weight 

84.32 

100.09 

147.63 

197.35 

Den.sity at 25^C, kg/m^ 

2963 

2707 

3700 

4437 

(Ib/ft®) 

(185) 

(169) 

(231) 

(277) 

Dissociation Constant, 
K, , atm 
■ 727 »c 

539 

6.4 X 10 "^ 

3.1 X 10~2 

4.6 X : 

927"C 

— 

1.8 

2.9 X 10 “ 

6.5 X : 

Technical-Grade 

Salt Cost, $/kg ($/lb) 

0.22 

0.07 

0.49 

0.29 


( 0 . 10 ) 

(0.03) 

( 0 . 22 ) 

(0.13) 
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The alkali and alkaline earth carbonates of Interest as phase'-change 
TES materials are readily available In bulk quantities as technical-grade 
salts, costing from $0. 07/kg to $2. 05/kg ($0. 03/lb to $0. 93/lb), as 
indicated in Tables 10 and 11. 

Table 12 contains a complete listing of the carbonate-containing salts 
initially considered as candidate salts, with the experimental or estimated 
values of properties relevant to their applicability as heat of fusion TES 
materials. These salts were selected for consideration because thier melting 
points lie within the general 343* to 871*C (650* to 1600*F) range, yet they 
possess reasonably high heats of fusion. Materials judged to have toxicity, 
cost, or availability problems were excluded from consideration. The 
materials considered include single salts and binary and ternary salt mixtures 
of both eutectic and non-eutectic compositions. 

The candidate carbonate compositions were screened at two levels. First, 
differential scanning calorimetry (DSC) was employed (when possible) to 
determine the melting and solidification temperatures and heat of fusion for 
the mixtures. Since the DSC's maximum operating temperature is 650*C (1200“F) , 
DSC was used only for steam TES applications. Then, based on heat of fusion, 
stability, and economic considerations, the most promising candidates were 
tested in laboratory-scale TES modules to determine the cyclic charge/discharge 
characteristics and short-term stabilities. 

Table 13 lists the salt mixtures selected for lab-scale testing, with 
their pertinent thermophysical properties. Primary emphasis for the salts 
selected for lab-scale investigations was given to alkali carbonate-based 
salt compositions with melting points in the solar-thermal range of 538“ to 
871*C (1000* to 1600“F), although for comparative purposes, fluorides and 
chlorides, when subjected to identical handling and fabrication procedures, 
were also considered. 


Perkin-Elmer DSC-2 differential scanning calorimeter. 
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Table 12. CANDIDATE CARBONATE-CONTAINING SALTS ILWING MELTING POINTS IN THE 

RANGE 398“ TO 898“C (750“ to 1650“F) 


Srtl t iMU 

Couij^auU Ion 

>loi t ' [ ~ ]h' X _ 

MoltlnH I’otut 

Huat of 

J/M X W 

FuHion 

Btu/Jh 

Coat 

"izji 

M /I’O 00 i 

Ml- l.'-.’o 

12-11- 15 

107 

74? 

2.77 

U9 

0,84 

0.38 

KjOO j-MsCO j 


68-12 

460 

8b0 



0.35 

0. 16 


(U-17 

47-51 

488 

010 

3.01 

168 

1.10 

0,54 

l.l iCO 

S UA7 

44-56 

4«6 

025 

4.Ui 

170 

0.05 

0.43 

U^l’O K;iCO) 

M 1-t7 

28- ? 2 

408 

028 

3.14 

1 15 

0.88 

0,40 

■ 2^'.0^ 


15-6 T 

505 

041 

1,42 

\u7 

0.00 

0.45 

KO l^-NiHCO pNa >S0,'j 


42,ii-14.‘). 17. 7 

544 

ion 

2.67 

m 

0.00 

0,04 

Kin-K^iH0/rNH^u‘0^ 

4?, \- 14-18. 7 

14 , 6-46 .0- 1 ‘J , 4 

5/0» 

1015 

2.74 

UB 

0.24 

0.11 

Li ^00 \“Na20O^-Ky,'O ^ 

20 -no- 20 

20-60-20 

5 50 

1022 

3.20 

IIH 

0.55 

0.25 

Lt^u’O pN.‘» u’O j-K >00 j 

22-lh-i>2 

22“ 16- 6 2 

550 

1022 

1.00 

120 

0.73 

0,33 

K01-K;00^«*Na^>00 j 

42.S-.. 1. 1-18.1 

28, 6 -5m. 0-17. 4 

558 

10 1/ 

2, 56 

110 

0.26 

0.12 

KjOiVpNaK“Na^>00 j 

42- 12- 2 n 

58,6-1 1,6-2 7, 8 

562 

1044 

3.02 

130 

0.35 

0.16 

HaOOi-NaOl-NinCO^ 

21-22-17 

16.1-11,2-52. 7 

562 

1044 


— 

0,37 

0.17 

K^OO^-NaK 

4^-5 > 

72.^1-27, 1 

168 

1055 

3.6 3 

156 

0.51 

0.23 

CaOO pOaF j-Ca (010 ^ 

2‘). 7-20. 1-50.2 

16.0-1 ‘1,0- 4 5.0 

57 ) 

106? 



... 


NaOl-NaK~Na uX> j 

42.5-20,5-17 

14. 2-11.8- 54,0 

575 

106 7 

4 , 04 

174 

0,13 

0.06 

OaOOpMF 

2‘). ?“70, 1 

62.0- 18,0 

»8U 

Ui;6 

-- 


2.14 

0,9 7 

K01*>Na a'O ^ 

50 , 2 -4‘) , 8 

4 1 . 5-58, 5 

588 

Ulo;< 

3.02 

nn 

0.07 

0.03 

LiF-l.i^OO) 

4‘)-5l 

25. 2“ 74. 8 

608 

1126 

7.16 

308 

2.03 

1,33 

NaOl-Na^SOfpNa|iOO \ 

51,a*-24. 1~24. 1 

1 1.6- 18.0-28.4 

61.' 

1114 

1.00 

120 

0,18 

0.08 

OaOl^-OaOOi 

70- 10 

72, t-27.8 

622 

U »2 


— 

0.00 

0.04 

KO l'-K^ CO j 

M.h- 18.4 

4 1> , - 5 1 . (> 

6 12 

1 1 ?0 

2.72 

11? 

0.24 

0.11 

Na01~Na^‘00 \ 

7t».'>-.2 1. 1 

<>4, 7“ l5. 1 

6 t.'f 

n?i 

4.05 

174 

0.04 

0.02 

KC1~K^h:o 1 

m- 15 

50-50 

6 16 

n;; 

2 . 70 

120 

0.24 

0.11 

OaCOpCa(OU);: 

15.8-64.2 

4 l.;l-5 7,0 

6 5 1 

120? 



0.09 

0.04 

Crj|Oj-Na^u\l j 

45-55 

54.0-46.0 

658 

1216 



1.81 

0.82 

OaOOpl^ljOO^ 

17-6 1 

44. 1-55. 7 

66.’ 

1224 




1.17 

0.53 

SrOl^-SrOO j 

74-26 

75. 1-24. 7 

671 

1240 



.... 


KF-K^OO ^ 

60-40 

18. 7-61. 1 

678 

1252 

3.1 2 

134 

0,84 

0,38 

llaOO pNa a*0 j 

17-61 

52.2-4 7.8 

686 

126? 



0,18 

0.08 

KF~KtOU ) 

7»6-54 

26.4-71.6 

688 

1270 

2 / 

110 

0.71 

0.32 

NaF“*Na jCO } 

18.7-61. 1 

20.0-80.0 

600 

1274 

3.72 

160 

0.20 

0.09 

Na;>00 pNa ,0 

61. 1- 18.7 

71.0-27.0 

605 

1283 

3.05 

170 


.... 

K;>00 pNa>CO \ 

18 .5-61.5 

44. ‘1-55. 1 

605 

1281 

2. 37 

102 

0.22 

oao 

K^>00 pNa >0O j 

44-56 

50.6-40.4 

710 

1110 

2,32 

100 

0,24 

0.11 

LIjOOj 

100 

100 

721 

nil 

6,07 

261 

2.05 

0.93 

UatUpHivOO 3 

51-4? 

54. 1-45,7 

714 

1151 



0.31 

0,14 

BaOO pNaOl 

47-51 

75.0-25,0 

714 

1353 



0.24 

0.11 

OaOO pNa 00 pNa ^ 80/, 

17,5-46.5-16 

14. 1-45.0-20. 7 

770 

1418 



0.11 

0,0.5 

OaOO pNa^HnVpNa,)SO,,j 

18- 18-24 

11.8- 15. 8-10.4 

705 

1461 



0,15 

0.07 

OrjOpK2<-'Vl 

20-80 

21.6-78.4 

80 1 

14 78 



1,04 

0.47 

Na^OO j’-Na^FO/, 

66.7-11. 1 

50.0-40. 1 

812 

1404 

2.26 

07 

0.18 

0.08 

HaCl^-llaCOj 

76-24 

86,8-11.2 

814 

1407 

— 


0,33 

0.15 

Haa\pNrttU 

24-76 

15.1-64.0 

814 

1407 

... 


0.13 

0.06 

Na^CO .j-NajBO/, 

62-18 

54.0-45. 1 

826 

1510 

2.21 

05 

0.20 

0,09 

Na^OOj 

100 

100 

858 

1564 

2,58 

111 

0.07 

0.03 

HaOl^-BaOO^ 

82.5-17.5 

81.1-16. 7 

860 

1580 


.. 

0.31 

0.14 

K^O.O \ 

100 

100 

808 

1648 

2.00 

86 

0.42 

0,19 
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System Components; Construction, Operating Procedures and Data 
Collection 

Construction 

Thirteen lab-scale containers, 12.7 cm (5 In.) high, with a 7.62 cm 
(3 in.) OD, and with a 1.27 cm (0.5 In.) OD single-pass internal HX tube, 
were fabricated according to the drawing shown in Figure 19. The outside wall 
of the cannister was Type 304 SS; the HX tube was Type 316 SS. The cannister's 
top was originally designed as a 1-ln. -thick translte disk bored to fit 
over the central HX tube; the thermocouple (T/C) sheaths were initially 
installed when 3.18 mm (1/8 in.) holes were drilled in the top and 
high-quality 4-mm OD, closed-one-end alumina tubes (Ventron Degussit) 
were inserted slowly into the molten salt. Later this method and design 
were found unsatisfactory for high-temperature (above 650®C) testing because 
transite deterioration, salt creepage, and alumina sheath thermal cracking 
occured. Therefore the translte cover was replaced with a Type 316 SS 
cover machined from 6.35 mm (1/4 in.) plate stock that had two holes for SS 
316 T/C sheaths and one vent hole. Swagelok fittings were butt-welded onto 
the cover to seal the stainless-steel sheaths. This arrangement performed 
well. After the salt mixture was loaded into each cannister, heli-arc 
welding secured this top to the cannister lip and OD of the HX pipe. A 
photograph of the final design is shown in Figure 20. Standard 1.27 cm 
(0.5 in.) Swagelok fittings connected the HX pipe to the TES system's air 
supply tubing. 

System Operating Procedures and Data Collection 

To facilitate lab-scale TES testing, a four-station test stand was 
constructed, with cooling air supplied through 1.27 cm (0.5 in.) diameter 
tubing manifolded to each station from a single source (Figure 21) . All 
lab-scale salt containers were inside two clam-shell resistance heaters 
(Thermal Corp., 825 watts, 1204“C, 2200°F) used to charge the units. The 
containers and heaters were wrapped with an 20.32 cm (8 in.) radial thick- 
ness of Johns-Manville Cerablanket high-temperature insulation (128 kg/m^, 

8 Ib/ft^ density), rated at a nominal thermal conductivity of 0.12 W/m-K 
(0.07 Btu/hr-f t-“F) at 540“C (1000°P) . The thermal barrier beneath the 
the unit consisted of 3.81 cm (1.5 in.) of firebricks, 12.7 cm (5 in.) of 
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Figure 20. FINAL DESIGN OF TES CANNISTER 
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insulation and 6.35 cm (2.5 in.) of firebrick as a base. A 17.8 cm (7 in.) 
thick layer of insulation was also used to minimize radiant heat losses from 
the top of the unit. A schematic of a TES module, comprised of the TES 
unit, heater, insulation, instrumentation and controls, is shown in Figure 22. 

Prior to welding of the top cover and installation into the module setup, 
each cannister was final filled with approximately 0.91 kg (2 lb) salt 
powder of the desired composition. Filling was accomplished by the following 
procedures. First, each container was filled with powdered salt, placed in 
a HOTPACK resistance box oven, and heated to 50‘’C (90°F) above the salt melting 
point. This step was repeated until approximately 80% of the container's 
volume was filled with molten salt. 

The units were operated under varying cycle periods (a cycle consisting 
of charge, holding time, and discharge). All salts studied in these units 
were cycled from mp + 50‘’C (90®F) to mp - 50°C (90°F) for a minimum of 
20 cycles, and they were maintained at these operating temperatures from 
228 to 1,032 hours. 


The salts were discharged with ambient air regulated by a Dwyer flowmeter 
[range 0 to 11.3 m^/hr (400 ft^/hr)]. The initial condition for a discharge 
run was that all T/C's, which were located at mid-salt level, would attain 
mp + 50°C (90°F). At the end of the holding period, power to the clam-shell 
heaters was cut, and air was discharged at a constant flow rate, which 
typically varied from 1.42 to 11.32 m^/hr (50 to 400 ftVhr). Mid-salt 
and air temperatures were continuously recorded on an Esterline Angus 24-point 
strip chart recorder throughout the cycle. A discharge run was ended when 
the slowest cooling salt T/C reached the salt's mp — 50°C (90“F) temperature. 
At this point, air flow was shut off, and the system was re-charged for the 
next run. The charging sub-cycle was controlled by a T/C at salt mid-level, 
between the cannister 's OD and the heater coils. The temperature controllers 
were set to bring this location to the holding temperature (mp + 50°C) in 
less than 2 hours, while initial heat-up rates were maintained between 3.5° 
and 10°C (6° and 18°F) per minute. 

After the end of each test, the containment and salt materials were 
examined to evaluate their stabilities. The salt composition was chemically 
analyzed and verified both before installation and after shutdown for CO 2 
content, cation composition (Li, K, Na, Ba, Ca, etc,), and relative levels 
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NO. THERMOCOUPLE 

cm 

in 


SALT 

1.12 

0.44 


SALT 

2.06 

0.81 


SALT 

3.02 

1.19 

(S) 

AIR OUTLET 

0 

0 




I7.8cm(7ln) of 

128 kg/m3(8lbs/ft3) 

(Stockad) 

RESISTANCE HEATER 
(Half Shall) V 


20,3cm (8 in) of 
128 kg/m®(8lb8/ff®) 

INSULATION 

(Radio%) 

20.3cm 


3.l7cm(l.2Sin)Hi 


12.7cm (Sin) of 
128 kg/m5(8lb8/ft®) 
( Stacked ) 


6.35cm (2.5in) High 


AIR FLOW, OUTLET 


INSULATION 

— MID- SALT 

RESISTANCE HEATER 
p(Half Shell) 


I7.8cm(7in) 


20.3cm(8in) of 
128 kg/m3(8lb8/ft®) 

INSULATION 
( Radially) 

20.3cm 


FIREBRICK 


INSULATION 


FIREBRICK 


3.18cm 1.25 in 


6.35cm (2.5 in) 


' — AIR FLOW, INLET 

Figuri' 22. SCHRMATIC OF LAB-SCALE MODULE 
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of Iron, nickel and chromium (elements Indicative of corrosive reaction 
between salt and stainless-steel ccnLainment vessel). Finally, each unit 
was sectioned at mid-salt level and visually examined to compare salt 
crystal formation patterns and to note overall compatibility. Selected 
HX tubes were sectioned lit the miniscus region of the salt, mounted, and 
metallurgically ejcamined to determine the nature and extent of corrosion. 

Salt TES Per f ormance Evaluation 

C riteria for Pe rf ormance Evaluation 

Numerous factors must be considered in any evaluation of TES system 
performance. Thus, a particular evaluation should be based on those factors 
most influential in TES system performance in a specific application (such 
as solar-thermal power system Integration) . The following parameters were 
used in comparatively evaluating the thermal performance of the salts in the 
lab-scale systems: 

• Average discharge heat flux during a discharge run 

• Variation of the discharge heat flux with time during a discharge run 

• Discharge salt solidification temperature range and characteristics 

• Cycling stability - possible variations in discharge behavlor 
resulting from salt/materials instabilities during thermal cycling 

• Operational stability — post-test examination of salts’ and 
containment materials’ conditions. 

Although several other parameters could be. added to this list, these 
form the basis for an assessment of a salt’s TES feasibility and for a rela- 
tive impact analysis of design modifications (such as thermal conductivity 
enhancement or alternate HX designs). 

As with the engineering-scale testing, heat fluxes were calculated using 
experimental data from Equation II. The heat-transfer area used to calculate 
this heat flux was defined as the portion of the central cooling tube in direct 
contact with the salt. Also, the air outlet temperature was measured at the 
center of the IIX tube near the salt level to yield heat fluxes that are not 
influenced sigtiificantly by radiant heat loss in the extended portions of the 
HX tube. 
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Evaluation of Salts for Steam-Ranklne Applications 343° to 399°C 
(650° to 750°F) 

As discussed in Task 1, the use of a TES subsystem to provide steam 
superheating can significantly increase the overall efficiency of a steam-Ranklne 
plant. For example, increasing the temperature of the steam entering the turbine 
from the TES boiler from 316° to A55°C (600° to 850°F) was estimated to result 
in an increase in system efficiency of from 20% to 30%. 

The conceptual TES unit for providing 455°C steam utilized four 
superheating stages, with salt melt temperatures in the range 349“ to 479°C 
(661° to 894°F) . Melting points in this general temperature regime can be 
attained with mixtures based on LI 2 CO 3 , Na 2 C 03 , and K 2 CO 3 . The lowest -melting 
carbonate eutectic known is the ternary composition 43.5 mole percent (mol %) 
Li 2 C 03 , 31.5 mol % Na 2 C 03 , 25 mol % K 2 CO 3 , As shown in its phase diagram 
in Appendix B ''which contains other available carbonate systems) , this eutectic 
has a melting point of 397°C (747°F) and thus is of interest for steam 
superheating applications. Also experiments in the U.S.S.R.^^ have indicated 
that additions of alkaline earth carbonates (5 mol % MgC 03 , CaC 03 , or BaC 03 ) 
to the ternary eutectic result in quaternary mixtures with solidus temperatures 
in the range 355° to 370°C (670° to 700°F) . Workers at Battelle Memorial 
Institute reported^ ^ that a 19.1 mol % addition of CaC 03 to the ternary eutectic 
resulted in a mixture having a solidus of 335°C (635°F) and a liquidus 
temperature of 481°C (898°F) . Such a mixture would allow heat storage and 
removal over a range of temperatures close to those identified earlier for 
the 455°C steam case. 

In an attempt to verify these reported melting point depressions, 
the thermal analysis technique of DSC was used to determine whether a 
lower melting (<397°C) carbonate will result from the addition of alkaline 
earth carbonates to the (Li-Na-K) 2 C 03 ternary. Separate additions of 5, 10 
and 15 mol % of CaC 03 and 5 mol % of SrC 03 were made to the ternary eutectic 
for this study. Upon heating at a rate of 2.5°C (4.5°F) per minute, the 
ternary eutectic began to melt at 397°C (747°F), in agreement with the 
published mp value, and melting was completed at 403°C (757°F). Upon 
cooling at the same rate, a small degree of supercooling was observed; 
solidification started at 387°C (729°F) and was completed at 384°C (723°F). 
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No supercooling effects have been observed with this ternary eutectic in the 
laboratory-scale TES unit. DSC determinations for the ternary and quaternary 
samples with alkaline earth additions are summarized in Table 14. An addition 
of 15 mol % CaC 03 to the (Li*-Na-K) 2 C 03 ternary lowered its mp from 397® to 
383®C (747® to 721®F), and an addition of 5 mol % SrC 03 brought the 
ternary's mp down to 389 ®C. 


Table 14. DSC RESULTS FOR THE (Li- Na-lOaCOs TERNARY AND TERNARY WITH 

ALKALINE EARTH CARBONATE ADDITIONS 




Melting 

Point 

AH^ 


Salt System 

Wt % 

"C 

°F 

J/kg 

2.67 X lO'" 

Btu/Xb 

1 . (Ll-Na-K) 2 C 03 
Ternary 

32-33-35 

396.2 

745.2 

114.7 

2. (Ternary) + CaC 03 
(5 mol’ % CaCOs) 

95-5 

392.8 

739.0 

1.95 X 10^ 

83.7 

3. (Ternary) + CaC 03 
(10 mol % CaCOs) 

90-10 

389.4 

732.9 

2.51 X 10*5 

108.1 

4. (Ternary) + CaC 03 
(15 mol % CaC 03 ) 

85-15 

383.0 

721.0 

1.90 X lO'^ 

81.8 

5 . (Ternary) + SrC 03 
(5 mol % SrC 03 ) 

92.8-7.2 

389.4 

732.9 

2.52 X 105 

108.2 


These results demonstrate the same general trends as those reported by 
the Soviets and Battelle Memorial Institute. While the 5, 10 and 15 mol 7 « 
CaC 03 additions exhibited distinct melting points, the Battelle results 
indicated that a 19.1 mol % CaC 03 fiddltlon resulted in a system having solidus 
and liquidus temperature 146°C (263°F) apart. Further investigation of the 
melting behavior of the (Li-Na-K) 2 CO 3 eutectic with alkaline earth carbonate 
additions appears to be needed. 

32 wt % LI 2 CO 3 - 33 wt % Na 2 C 03 - 35 wt % K 2 CO 3 (Module L3-1) . 

To further investigate the ternary's applicability to steam utility TES 
subsystems, the thermal behavior of the salt was studied in a lab-scale 
cannlster with the original transits top design, which was found to perform 
well for this low-temperature salt mixture. In this initial design, salt 
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temperatures at mid-depth were monitored during discharging with room-temperature 
air, using Type K (chromel/alumel) T/C's situated at three different radial 
locations of r ^ i.i, 2 .I and 3.0 cm (0.44, 0.81, and 1.19 in.). These 
radii correspond to reduced salt thickness values (r/r^-1) of 0.76, 2.24, 
and 3.76, which are used in the Megerlln-type heat transfer analysis.^** 

Based on the fact that the volume of salt contained between the HX tube and 
radius (r) in the container is proportional to r^ - r^^ (where r^ = 6.35 mm 
or 0.25 in. HX radius, r = T/C radial distance), these T/C locations correspond 
to solldlflcatlons of 7%, 33% and 74% of the total salt volurae when the 
solid/liquid interface reaches each respective location. 

Module L3-1 was loaded with 611 grams (1.35 lb) of the ternary eutectic 
salt. Characteristic discharge curves were determined as a function of 
discharge air flow rate over the range 1.42 to 11.32 m^/hr (50 to 400 ft^/hr). 

Figure 23 contains a plot of the discharge time required for the salt to 
solidify versus air flow rate. At air flow rates greater than 4.25 m^/hr 
(150 f t^/hr) , the rate of discharge appears to be controlled mainly by heat 
transport through the solidified salt. A flow rate of about 5.1 m^/hr (180 ft^/hr), 
which corresponds to a linear flow velocity of 15.24 m/ second (50 ft/ second) 
through the HX tube, was chosen as the standard, optimum flow condition. This 
flow rate represents a Reynold's number of 2.1 X 10**, thereby satisfying 
the turbulence, criterion. A similar determination of the optimum discharge 
flow rate was carried out for all salts studied in the lab-scale TES cannlsters. 

Results obtained during a typical run (No. L3-1-3) of the salt temperature 
discharge profile for an air flow rate of 5.1 m^/hr (180 ft^/hr) are plotted 
in Figure 24. From the time-temperature response of the salt T/C's, it appeared 
that the salt was solidifying congruently over the temperature range 410“ to 
400“C (770° to 752“F) taking 75 minutes for 74% of the total salt volume to 
solidify. Figure 25 shows the Instantaneous heat flux curve calculated for 
this salt as a function of discharge time for an air flow rate of 5.1 m^/hr 
(180 ft^/hr). From this, the average discharge heat flux of 6.95 X 10^ W/m^ 

(2,206 Btu/hr-ft^) was graphically determined, as shown in Table 15 summarizing 
the experimentally determined data. 
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Figure 23. TYPICAL DISCHARGE PERFORMANCE OF LAB-SCALE MODULE L3-1 
CONTAINING TERNARY (Li-Na-K) -CO- EUTECTIC 
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After 10 cycles and 504 hours at operational temperatures between 450“C 
(842®P)(mp + 50®C) and 350®C (662“F)(mp — 50*C), the unit was shut down and used 
as the TES medium for an analysis of a porous aluminum TCE material that will 
be discussed in detail In Section 2.3 of this report. The data generated 
during these 504 hours were utilized to establish a performance baseline for 
further comparative testing. For this Initial set of tests the unit operated 
with excellent stability, no loss of salt due to creepage, and no equipment 
failure. 

Evaluation of Salts for Solar-Thermal Applications 538® to 871°C 

(1000® to 1600®F) 

Advanced concepts for the conversion of solar-thermal energy to electrical 
power Include large central receiver systems employing a closed Brayton 
power cycle, as well as small dispersed power systems with Stirling heat engines. 
TES subsystems are desirable In these applications to provide energy to the 
turbine or heat engine during periods of low solar insolation - early morning, 
late afternoon, and periods of cloud cover. Latent heat storage Is particularly 
attractive because of Its potential for attaining high-energy storage densities 
and its ability to receive and release heat over a narrow temperature range. 

The salts investigated in this subtask for their applicability to latent 
heat storage were in the temperature range 538® to 871®C (1000® to 1600®F) , 
which represents the expected regime of storage temperature requirements for 
storage-coupled advanced solar systems.^^’^®’^^ These high design temperatures 
are being used to ultimately obtain higher thermal conversion efficiencies. 

The following account details the operational life of the first three 
high-temperature salts selected for lab-scale evaluation. It was during 
these Initial tests that some of the significant problems associated with 
molten salt systems in this high-temperature regime were encountered and 
solutions identified. The first salts tested were 49.4 wt % Na 2 C 03 - 50.6 
wt % K2CO3; 1.21 wt % LI2CO3 - 49.95 wt % Na2C03 - 48.84 wt % K2CO3; and 
pure LI2CO3. 

49.4 wt % Na 2 C 03 - 50.6 wt % K 2 CO 3 (Module Ll-1) . This 
Na 2 C 03 -K 2 C 03 mixture was expected to correspond to a congruently melting 
salt with a melting point of 710®C (1310®F), tlie lowest-melting composition 
in this binary system. The phase diagram for this, and other available 
carbonate systems, are included in Appendix B. 
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A cannlster was loaded with 880 grains (1.94 lb) salt of composition 
49.4 wt % Na 2 C 03 - 50.6 wt % K 2 CO 3 and assembled with the original cannistei 
design of transits cover and high quality alumina T/C sheaths as previously 
described. The molten salt level in these lab-scale units was approximately 
1.27 cm (0.5 in.) below the top of the vessel; the salt volume in a cannlster 
was approximately 426 cm^ (26 In.^). 

As for the ternary, and for all salts studied in this subtask, characteristic 
discharge curves were Initially obtained as a function of discharge air flow 
rate over the range 1.42 to 11.32 m^/hr (50 to 400 ft^/hr) to determ' ,c the 
optimum flow condition for that unit, so that the rate of discharge (salt 
solidification) is controlled mainly by the heat flux through the solidifying 
salt. The plot of salt solidification time versus air flow rate appears in 
Figure 26, from which 180 ft^/hr was chosen as the standard flow rate for 
this system. 

Discharge results obtained during a typical run (No. Ll-1-2) for an air 
flow rate of 5.1 m^/hr (180 ft^/hr) are plotted in Figure 27. This time- 
temperature response of salt T/C's confirmed that the salt was indeed 
solidifying congruently over the temperature range 708® to 701®C (1306® to 
1294®F) . No significant thermal gradients were observed through the still 
molten salt regions during solidification, indicating strong convective mixing 
effects consistent with previous TES results. Figure 28 shows the instantaneous 
heat flux curve calculated for this salt as a function of discharge time for 
an air flow rate of 5.1 m^/hr (180 ft^/hr) . 

After 13 charge/discharge cycles and 336 hours of operation of Module Ll-1 
with the Na 2 C 03 -K 2 C 03 salt, the heaters became nonfunctional. Upon disassembly 
of the unit, it was discovered that a significant amount of salt had crept up 
the inner container wall and outer HX tube surfaces, reacting with the 
insulation and corroding the metal furnace 1 Js. The transite cover initially 
used over the container had also deteriorated as a result of the severe 
temperatures and chemical environment. The color of the solidified salt had 
changed from white to blue, possibly from corrosion of the stainless-steel 
container. Because of these experimental difficulties, the following 
cannlster design modifications were made: 
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Figure 28. EXPERIMENTAL HEAT FLUX CURVE FOR MODULE Ll-1, 49.4 wt % Na2C03 — 50.6 

AT 5.1 m3 /hr (180 ft^/hr) AIR FLOW RATE 
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1. The translte cover was replaced with a machined stainless-steel cover, 
containing three holes for T/C sheaths and one vent hole. This cover 
was welded onto the container wall and HX tube, greatly reducing the 
area that allowed salt to creep out the container, yet still utilizing 
the alumina T/C sheaths. 

2. Lost Na2C03-K2C03 salt was replenished, and 1 mol % LI2CO3 was added 
to this same cannlster to reduce corrosion of stainless steel by 
possibly forming a LlCr02 protective layer on container and HX walls. 

This unit was then run as Ll-2. 

1.21 wt % LiaCOa - 49.95 wt % NaaCOa - 48.84 wt % K2CO3 (Module 
Ll-2 . After the above modifications were completed, the cannlster contained 
888 grams (1.96 lb) of the salt mixture and was operated as Module Ll-2. The 
discharge and heat flux curves obtained at a flow rate of 5.1 m^/hr (180 ft^/hr) 
are shown in Figures 29 and 30 i*espectlvely . From the tima-temperature response 
of the salt T/C*s, it appeai'ed that the doping of the Na-K system with Li2C03 
slightly altered the melting behavior of the system by causing it to solidify 
incongruently over the temperature range 710° to 695°C (1310° to 1283°F). 

Figure 31 shows the variation of heat flux with discharge time for four cycles 
at different points during the lifetime of this unit. The genex'al shape of 
the curves remained unchanged while defining a narrow lifetime discharge 
performance band, indicating no significant loss of thermal storage capacity 
or hear transfer capability during more than 500 hours of operation. 

Problems were encountered with thermal cracking of the recrystallized 
alumina T/C protection sheaths after 140 hours of operation. Failure of the 
alumina sheaths allowed salt contact with the 0.51 mm (0.02 in.) thick-walled 
Type 304 SS clad T/C's and led to eventual T/C failure. Because ceramic sheath 
failure was a recurring problem, the T/C design was modified to incorporate 
three 0.76 mm (0.03 in.) thick-walled Type 316 SS sheaths that were directly 
welded to the stainless-steel cover. 

This unit was restarted and finally terminated after 23 cycles and 528 
hours in operation, with no loss of salt or change in chemical composition 
of the mixture. 
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Figure 29. TYPICAL DISCHARGE PERFOIMANCE OF (Na-K)_CO, SYSTEM WITH Li^CO- ADDED 




LAB MODULE 



Figure 30. EXPERIMENTAL HEAT FLUX CURVE FOR (Na-K)-CO- SYSTEM WITH Li„CO^ ADDED 
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Figure 31. LIFETIME TISCBARGE PERFOBMANCE OF (Na-K)-CO_ SYSTEM WITH Li_CO ADDED 
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L12C03 (Modules L2-1 and L4-2). Originally Module L2-1 contained 
720 grains (1.6 lb) of LI2CO3 and operated with the transite top design during 
the same period as the Li-doped (Na-K)2C03 system. 

Typical discharge and heat flux curves are shown in Figures 32 and 33, 
respectively. Solidification Oi the L12C03 was observed to occur within 
+ 3°C (5.4°F) of its quoted mp of 723‘’C (1333‘’F). The salt at the 
outer T/C location solidified before that at the middle T/C because of heat loss 
through the outer container wall. 

This system encountered difficulties with alumina T/C sheath failures 
simulataaeously with the Li-doped (Na-K)2C03 system. After 13 cycles and 
408 hours the unit was shut down for post-test examination. Visual examination 
revealed salt leakage through a small pinhole in the bottom weld and a similar 
defect in the container wall at about mid-salt level. Post-teit analysis of 
the salt revealed an Increase in iron content from 0.03 to 0.36 wt % and an 
increase in nickel from >0.01 to 0.03 wt %. 

Based on the successful results obtained when the design of the Li-doped 
(Na-K)2C03 system was modified, a new final can design was established. This 
final design Incorporated two 10.22 mm (0.4 in.) thick-walled open-one-end 
Type 316 SS sheaths and 3.18 mm (0.125 in.) OD T/C's with 1.02 mm (0.04 in.) 
thick Type 304 SS cladding. It avoided direct welding of the T/C sheaths by 
using butt-weld tube fittings (Swagelok) , with a welded vent tube that could 
also serve to provide an environment of CO2 above the salt (as in Figure 20) . 

This design proved generally satisfactory for subsequent laboratory-scale 
testing. 

In Module L4-2, 754 grams (1.7 lb) of Li2C03 was retested in the modified 
cannlster described above with a cover-gas environment of CO2 over the melt. 

From the salt discharge plot of Figure 34 one can see that this system 
solidified between 734° and 730°C (1353° to 1346°F). Also, from the instantaneous 
heat flux profile of Figure 35, the average heat flux to the air was calculated 
to be 76,574 W/m^ (24,290 Btu/hr-ft^) . This was the largest heat transfer 
exhibited by any salt, regardless of temperature regime. The consistency of 
this unit's performance becomes evident on inspection of Figure 36, which plots 
the discharge heat flux from four different cycles in the unit's lifetime. 
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Figure 32. TYPICAL DISCHARGE PERFOBMANCE OF LiaCOa SYSTEM (MODULE L2-1: AIR COVER GAS) 

AT 5.1 m3 /hr (180 ft^/hr) AIR FLOW RATE 





Figure 33. EXPERIMENTAL HEAT FLUX CURVE FOR LiaCOs SYSTEM (MODULE L2-1; AIR COVER GAS) 

AT 5.1 m^/hr (180 ft^/hr) AIR FLOW RATE 
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Figure 36. LIFETIME DISCHARGE PERFORilANCE OF Li-CO- SYSTEM WITH CO, COVER GAS 
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These profiles not only define a narrow discharge performance band but also 
seem to exhibit a performance stabilization period during the first seven 
cycles. 

This unit was maintained at operating temperatures for 672 hours and 
22 charge/discharge cycles, at which time it was shut down and sectioned for 
post-test evaluations. The solidified mid-salt cross section is shown in 
Figure 37. The size and morphology of solid salt crystals in this and other 
sectioned modules could be correlated with the manner in which the final 
cooldown was performed. In this particular instance, the salt’s last discharge 
run was a free-cool run, in which the salt was allowed to cool from 50“C (90“F) 
above its mp to ambient with no nir-cooling discharge. From this 
photograph, one may observe the somewhat fine radially-oriented columnar 
crystals that formed near the outer edge. This growth pattern is the result 
of fairly rapid outward cooling through radiant heat loss to the ambient. 

From this outer edge inward, the crystals become larger and randomly oriented, 
a result of slower cooling. 

The new cannister design was successful in that no salt creepage was 
observed throughout the cannister' s operational life. From Table 16 one can 
see that this unit's post-test analysis revealed no increase in iron or nickel. 
Also, the CO2 cover gas aided in reducing the corrosion rate compared with an 
air cover gas, most likely by minimizing Li2C03 decomposition to Li20 and C02. 


Table 16. PRE- 

AND POST-TEST CHEMICAL 

ANALYSIS 

OF LX2CO3 

SYSTEM 


Chemical Analysis, Wt 

% Li 

col 

Fe 

,Ni 

Cr 

Pre-Test 

19.08 

79.26 

0.03 

<0.01 

<0.01 

Post-Test (Air) 

18.89 

80.23 

0.36 

0.03 

0.01 

Post-Test (CO2) 

19.39 

80.53 

0.02 

<0.01 

<0.01 


The following synoposes highlight the most important experimental results 
of the remaining six salts selected for possible solar-thermal applications. 
All utilized the final welded containment design. 
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Figure 37. TRANSVERSE CROSS SECTION THROUGH 
MODULE. L4-2 Li2C03 SYSTEM AFTER 672 HOURS AND 
22 CYCLES 
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Na2C03 (L2-2) . Although Na2C03 pOv«^sesses a lower heat of fusion 
(2.6 X lO'" J/kg, 113 Btu/lb) than LiaCOa (6.1 X 10® J/kg, 261 Btu/lb), this 
salt is attractive because its mp of 858“C (1576°F) approaches the upper 
end of the solar temperature range being studied, and it may become 
cost effective ($0. 07/kg, $0. 03/lb) in a large-scale solar-power system for 
which the mass of TES medium is not restricted. 

With this in mind, a lab-scale unit containing 827 grams (1.82 lb) of 
NaaCOa was tested in a welded container. This discharge profile of Figure 38 
shows that the salt nominally solidified at 865°C (1589°F), while maintaining 
an average discharge heat flux of 45,648 W/m^ (14,480 Btu/hr-f t^-) > as calculated 
from che heat flux graph of Figure 39. This curve shows the instantaneous 
heat flux for Na2C03 to be decreasing at a rate of 788 W/m^ (250 Btu/hr-ft^) 
per minute, as compared to 113003 's decrease of 451 W/m^ (143 Btu/hr-ft^) per 
minute. The lifetime performance of the unit exhibited relative stability 
throughout Run 17, with a decrease by Run 22, Post-test examination revealed 
slight leaks around the T/C fittings, which resulted in a loss of approximately 
15 grams (0.03 lb) of salt. 

The free-cool shut down occurred after 21 cycles and 288 hours of operation. 
Chemical analysis of the terminated salt indicated no change in nominal composition 
and 3 slight increase in iron (from 0.03 to 0.06 wt %) , a result of hot corrosion, 
as summarized in Table 17. When the. cannister was sectioned and photographed 
(Figure 40), micro-cracks and voids ranging from 0.025 to 6.35 mm (0.001 to 
0.25 in.) were found in the salts along very evident radially-oriented salt 
solidification patterns. These patterns are very different from the typical 
free-cool crystallization patterns found in the 113003 unit. If these voids 
and cracks were also formed in the solid phase under the "forced-cooling" of 
the standard operating procedures, they would act as thermal resistances and 
reduce the effective heat flux. If such defects were filled with air or CO2 
during discharge, they would have a thermal conducltivlty of 0.07 W/m-K 
(''^0.04 Btu/hr-ft^-°F) , compared to 1.73 W/m-K (^4 Btu/hr-f t-°F) for solid Na2003. 
These apparent factors could ultimately lower Na2003's effective thermal 
conductivity through the solid salt and partially explain the much lower heat 
fluxes observed for Na2C03, compared to 113003. 
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Figure 38. TYPICAL DISCHARGE PERI’OBMANCE OF Na2C03 SYSTEM 
AT 5.1 m3 /hr (180 ft^/hr) AIR FLOW RATE 
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Figure 40. TRANSVERSE CROSS SECTION THROUGH MODULE L2 
NajCOs SYSTEM AFTER 288 BOORS ASD 21 ctSs 
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Table 17. PRE- AND POST-TEST CHEMICAL ANALYSIS OF Na 2 COj SYSTEM 


Chemical Analysis, Wt % 

Ka_^ 


Fe 

_N.1_ 

^ Cr 

Pre-Test 

44 . 29 

56.40 

0.03 

< 0.01 

< 0.01 

Post-Test 

43.07 

57.39 

0.06 

< 0.01 

0.01 


Metallurgical analysis of the Type 316 SS HX tube at the meniscus of the 
salt-sioe surface revealed tlie formation of an outer, porous, hot corrosion 
layer and sub-scale carburization of both the bulk grains and grain boundaries 
as a result of salt interactions. Air-side oxidation of the HX was also observed. 
Additimially, precipitated particles of what appeared to be brittle inter- 
metallic sigma phase were observed throughout the bulk stainless-steel matrix. 
Sigma phase formation, a result of the prolonged high-temperature exposure and 
not of salt corrosion , is deleterious to the strength and related mechanical 
properties of stainless steel. Figure 41 shows a micrograph of this Type 316 SS 
which was exposed to NapCOs for 288 hours at 807° to 907°C (1485° to 1665°F) . 

Of particular note are the 150 to 200 pm ( 6 to 8 mil) thick porous, hot corrosion 
zone at the surface and the precipitated sigma phase particles within grains and 
grain boundaries. These clearly indicate that selection and testing of alternate 
containment materials are required for these high- temperature solar-thermal 
applications, not only on the basis of salt compatibility, but also on the basis 
of possible interactions with actual working fluids (air. He, molten sodium). 

K 2 CO 3 (L2-3) . Potassium carbonate, with a melting point of 898°C 
(1648°F), is the highest melting carbonate of interest in the upper temperature 
range for the solar-thermal applications. Therefore, to completely and thoroughly 
span this entire temperature range, this salt was selected for testing. 

This unit also utilized a CO 2 cover gas environment at the salt surface. 

Its discharge profile and its rapidly decreasing heat flux curve (Figures 42 and 
43, respectively) are highly suspect, though, because they represent only two 
cycle runs prior to system failure. A high discharge solidification temperature 
of 912°C (1674°F), and calculated average heat flux of 72,507 W/m^ (23,000 
Btu/hr-ft^) were observed, but additional confirmatory data would be desirable. 
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Figure 41. TYl’E 316 SS CONTAINMENT AFTER 288 HOURS EXPOSURE TO NajCOa AND 
21 CYCLES AT 807° TO 907°C (1485° to 1665°F) UNDER CO2 COVER GAS 

(FeCl3-HCl-H20 Etchant) 
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It appears that this salt defined the maximum useful operating temperature 
of AISl 304 and 316 for TES applications. After two cycles and 96 hours at 
948°C (1738°F), the HX tube failed, allowing no air to flow through it on 
discharge. As a result of tills, the unit was terminated prematurely and 
sectioned lengthwise along the cooling tube. The salt surface and embedded 
fragments of HX tubing are shown in Figure 44. Examination of this section 
revealed that the tube had become extremely oxidized and embrittled from 
the contlnous high-temperature exposure. This eventually led to a complete loss 
of structural integrity, resulting in a HX breach that allowed salt to enter 
and fill the HX tube. 

Post-test chemical analyses revealed that the detectable amount of hot 
corrosion products containing iron, nickel, and chromium were deposited in the 
bottom on the cannister. Lower concentrations of these elements were found at 
the top (Table 18) . 

Table 18, PRE- AND P08T-TEST CHEMICAL ANALYSIS OF K 2 CO 3 SYSTEM 


Chemical Analysis, Wt % 

_K 

CO 3 

Fe 

Ni 

Cr 

Pre-Test 

56.15 

43.05 

0.02 

<0.01 

<0.01 

Post-Test (Bottom) 

54.91 

42.76 

2.27 

0.03 

0.24 

Post-Test (Top) 

55.73 

44.28 

0.03 

<0.01 

<0.01 


47.8 wt % Na 2 C 03 - 52.2 wt % BaCOs (L3-3) . Very little thermo- 
physical property data were available for this system, other than that showing 
there was a congruently melting eutectic (reportedly at 686°C, 1267 °F)^-^ at 
this composition. As a prospective TES medium, this, and the next Li 2 C 03 -CaC 03 
system (mp '^^662'’G, 1224'’F), represent the only carbonates with melting 
points in the 650°C (1202 ®F) region. 

Since this system was attractive from an ecoiwmic standpoint ($0. 18/kg, 

$0. 08/lb), and little past work was conducted with it, 966 grams (2.1 lb) 
of this Na 2 C 03 ~BaC 03 mixture was cycle-tested for extended operations. It is 
important to note from Figure 45 that the thermal arrest in the salt's temperature 
versus time curve occurred at 717° to 712°C (1323° to 1314°F) instead of the 
published value of 686 °C. The consistency with which this solidification occurred 
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Figure 44. LONGITUDINAL CROSS SECTION THROUGH 
MODULE L2-3 K 2 CO 3 SYSTEM THAT FAILED AFTER 96 
HOURS AT 841° TO 941°C (1546° to 1726°F) AND 2 CYCLES 
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Figure 45. TYPICAL DISCHARGE PERFOSMANCE OF Ha 2 C 03 — BaCOg SYSTEM 
AT 5.1 (180 ft^/hr) AIR FLOW RATE 
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indicates that this Na2C03-BaC03 eutectic way be used for phase change TES 
applications in the higher 700® to 735®C (1290® to 1340°F) range. An experimental 
heat flux of 41,676 W/w^ (13,220 Btu/hr-ft^) was graphically calculated from 
Run 21 (Figure 46) . This plot also outlines the relatively stable discharge 
performance exhibited by this salt throughtout 27 cycles of its operating life. 

After 36 cycles and 984 hours of operation, this unit was shut down by 
forced cooling with air at a typical working flow rate (5.1 m^/hr, 180 ft^/hr) 
down to ambient. The sectioned photograph of Figore 47 illustrates the crystal 
size formation patterns characteristic of this type of discharge. As is expected, 
near the cooling tube the crystals are fine, with columnar shapes, exhibiting 
a distinct radial orientation with respect to the cooling tube. This type of 
growth pattern is the result of very rapid cooling. Midway between the cooling 
tube and outer wall, the crystals are larger and randomly oriented because of 
slovvrer cooling. Near the outer edge, the crystals are somewhat finer again but 
have no preferred orientation, indicating a higher cooling rate than that in 
the middle, but much slower than that near the HX tube. These observations are 
consistent with the temperature profile through the salt as a function of cooling 
time, as shown in Figure 45. 

Post-test analysis revealed no salt decomposition or significant containment 
corrosion. 

55.7 wt % Li2C03 - 44.3 wt % CaCO^ (L4-1) . With its reported 
rap of 662°C (1224°F)>^^ this was the lowest melting carbonate mixture tested 
for solar-thermal applications, A quantity of 814 grams (1.79 lb) of this 
mixture was loaded and run in a welded container, with extremely stable results. 
From Figure 48, it was observed that no super-cooling occurred and that the 
Ll2C03-CaC03 salt consistently solidified within i5®C (9°F) of its reported 
mp of 662®C (1224°F). 

Although the module contained 56 wt % 113003, its average heat flux to 
the air of 40,415 W/m^ (12,820 Btu/hr-ft^) is one of the lowest tested. This 
value was calculated from Run 25 (Figure 49, which also displays a very narrow 
■’ Ifetlme performance band for Runs 2, 12, 22, and 25.) 


LAB MODULE^ L3-3, — 69,344 

SALT= NogCOs-BoCOa 47.8 -52.2 wt% 



Figure 46. LIFETIME DISCHARGE PERFOSMAMCE OF Na2C03 — BaC03 SYSTEM OVER 984 HOLTiS 






CUTTING BL40E VOID, NOT 
SOLIDIFICATION VOID 


Figure 47. 


TRANSVERSE CROSS SECTION THROUGH MODULE L3-3 47.8 Z N«,CO, 
52.2 wt A BaC 03 AFTER 984 HOURS AND 36 CYCLES 
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Figure 48. TYPICAL DISCHARGE PERFOSMAircE OF L± 2 C 03 - CaCOo SYSTEM 
AT 5.1 m^/hr (180 ft^/hr) AIR FLOW PvA.TE 
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Figure 49. LIFETIME AVERAGE INSTAiTTASEOUS HEAT TUIK CLTvVE FOR Li CO -CaCX). SYSTEM 






S(uoo no hoal oi lusUm valuoa worn avaUalvlo In litoraturo for oither tho 


Na-'CXV^-BiUU)j or Iho Iil.'C'Oj-claCO^^ synloin, ost linalod vnluoK wore ol^talnod by 
oompariuy, Mu^ Inti'^rnt od (|uantUy of lioat romovod from tiu'so salts with that 
romovod i rom hi and Na,>CH) ^, wlu^so boats of fusion aro known. Tho rosults 
apiH'arod provlously In Tablo l b (paroulhol I cnl ly) . 

Moro rtvl tablo hoiM of ins ion and thormal oonduoLivily monsuromonts ari> 
roquirod to aoouratoly anaivno tho tl)ormal porformanco and altraoti.vonoss ol" 
tho and laf'CU) ^ systoms as PCM*s ior furthor TKS dovolopmont 

work ♦ 

Hits l.{.>Cclj-BahOj syslom was lorminatod in a froo-oool aflor 312 liours and 
2 b oyolt\s, l\\st--ti'st oxamlnaltcms rovoalod no sail doooinpos 1 1 ion and no Inoroaso 
(n conoonlral ion ol molallio oonlainiuonl oloinonts from oorrosion. 

sorviHl a tluai purposo durinp, its iab-soaio lostinp.. Although invostigatod as 
a cNUulldato TKS modium in Iho solar-lhormal 338^' to 871”C1 (1000^^ to IbOO'MO 
lompi'ratnro rango, its moro Important rami f icalions woro as a TCK oivncopL. 

A dotailotl aooounl i^f this unites oxporimontal rosults is prosonU'd in Si'otion 
2.3, tosllng of Ihormal oonduotivity onhanoomonL oonoojUs. 

Non-"tkn3unia ti\ Halls 

ll is a rocognlniHi fact that no si agio salt, or family of salts, possossos 
siiporior proporlios in all Important catogorios (siioh as low eost, liigli hoat’-of-' 
fusion, oaso of handling, low corrosivity, and high‘~tom]K'raturo stability)-, It 
is obvious tliat trado’^offs in thoso pnramotors aro nocossary whoa seloctiug a 
sal t for a givon appi ication and its part icular duty roqui rinonts . Thus it i s 
jvrobablo that carbouatt' salts may someday bo an optimum solt'cl ion for one 
application, with a fluorido, chloride, or other salt being a better choice 
for anot:1\er ap]vlication liaving different system performance and cost limitations. 

In keeping with this logic, and to estalilish a baseline of comparative 
thoroughness, two non-carbona to mixtures were tested in the lab-*scale eannlsters. 
Tlu\se salts were treated with tlie same handling, loading, and operating procedures 
that provetl successful for the carbonates. A closed cannister system was to be 
used, with an inert environment above the salt, as required. 
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GaCl;) (L 3--4) . This chloride system was chosen because its 772°C 
(1422°F) inp was compatible to the temperature requirements of the study, 
and its medium heat of fusion (2.56 X 10^ >T/kg, 110 Btu/lb) and low unit cost 
($0. 09/kg, $0. 04/lb) made it ai\ attractive TES candidate. 

Because of GaClo's hygroscopicity, the mixture was handled in an air dry- 
room and dehydrated overnight in a Uotpack oven purged with dry CO?. . Gravimetric 
analyses indicated nearly complete removal of the water of hydration. In keeping 
with the standard procedures, 727 grams (1.6 lb) of the salt were then loaded 
into a cannister and melted In the Uotpack oven in two fillings. In less than 
6 hours (at 822°C, 1512°F) in this oven, a large amount of salt crept out of the 
cannister. Upon iiimiediate removal and visual inspection of the cannister, it was 
found that excessive hot corrosion had scaled off a 180 pm ("^7 mils) layer 
rrom the CD of the. SS Type 304 cannister. 

The unit was finally filled to 2.54 cm (1 In.) from the top of the can 
with molten salt and maintaiixed at 725°G (1337“F) for 20 hours (overnight) , prior 
to welding of the top and closing of the system. Although this was often done with 
carbonates, the chloride cannister was found to have lost its operational integrity 
entirely, with a corrosion scaling loss of 635 pm (25 mils) from the cannister 
OD, and an eroded hole in the SS 316 IIX tube, as pictured in Figure 50. 

Unlike carbonates, the failure of this salt during loading demonstrated 
that chlorides require a dry, monitored Og-free inert exxvironraent throughout 
the entire handling, loading, and operating procedures to prevent extensive, 
corrosion of metallic containment. 

32.5 wt % N aF - 67.5 wt % KF (L4-3). This NaF-KF system is an 
attractive TES candidate because of its high heat of fusion (5.86 X 10^ J/kg, 

252 Btu/lb) and its melting point (721°C, 1330°F) . Also, it was selected as 
the most cost-effective fluoride for TES applications in the 700° to 725°C 
(1292° to 1337°F) temperature I'ange in a Pennwalt Corporation report by 
J. L. Eichelberger. 

As with the chloride salt, this salt \<ra.s to be operated eventually in a 
closed system under a high-purity argon environment. After two fillings and 
4 hours at 'v770°C (1418°F) , significant amounts of salt crept out of the containers 
and corrosively scaled off a 254 pm (vlO mil) layer of stainless steel from 
the outer surface of the cannister (Figure 51) . 
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Figure 50, MODULE L3-4 AFTER FILLING WITH CaCl 2 
AND HOLDING AT 725°C (1337°F) FOR 20 HOURS 
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Figure 51. MODULE lA-3 AFTER FILLING WITH NaF-KF 
AND HOLDING AT 770°C (1418°F) FOR 6 HOURS 
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AluTuinuTu powder was to be a **getter*' at this point, but too much of the 
camiister had corroded after this first 4 hours at operational temperatures 
to consider it functional. Once again, the necessity for a dry, Op-free 
environment above fluorides and chlorides to load and weld the containers became 
visibly apparent. 


If one considers these facts, it becomes clear that the application of either 
of these two non-carbonate TES mixtures on a large industrial scale would entail 
the use of sophisticated atmosphere control during loading, operation, and 
maintenance. The added expense for salt handling alone could significantly 
decrease the cost effectiveness of either the energy-dense fluorides or inexpensive 
chloride salts. These points, together with the uncertainty of personnel safety 
in the event of a system leak, become important institutional considerations of 
these salts’ ultimate TES applicability. 


Experimental Verification of Megerlln Analysis 


An accurate heat transfer model describing the melting and solidification 
of PCM salts is necessary if one is to adequately interpret the results from 
experimental TES modules and to meaningfully extrapolate these results to large- 
scale systems that use actual fluids at elevated temperature and pressure. The 
analytical model developed at the Institute of Gas Technology (IGT) during t?v.^ 
original contract was based on Megerlin’s analysis of the solidif cation of a 
cylindrical melt via an internal HX.^^ The general mathematical approach 
adopted was obtained by modifying some existing solutions for phase-change 
thermal behavior. One of the more appropriate solutions was that originally 
developed by Megerlin.^*^ Using a cylindrical geometry, a finite heat-transfer 
coefficient, and the liquid (salt) at its melting point, the time required to 
solidify a given thickness of salt can be expressed in a relationship containing 
a number of dimensionless parameters; 


N, 


Fo 


K 


rsVinr*(Ng^lnr* + 2) dr* 


(28) 


(Ng^lnr* + 1)2 + 2 ^ 


Bi 


ph 


Inr* (Ng^lnr* + 2) - (Nj^^lnr* + 1) 
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whei'e ” 

r . (the Fourier number) gives the dimensionless time required to solidify 
.( given thickness of salt, 

Npo “ (29) 

Ngi (Biot number) gives the ratio of the surface conductance to interior 
solid conductivity, 


“ ha/kg or Ua/k^ (30) 

Nph (phase change number) gives the ratio of the heat of fusion to the 
sensible heat between the melting point and a reference temperature 
(temperature of cooling/heating media), 


Pt 

P" ' °p ~ Pp (31) 

The term R is a reduced thickness; the ratio of the thickness of the annular 
layer of solid salt to the radius of the HX pipe. 

a = characteristic thickness, m (ft), (half thickness of a 
slab, radius of a cylinder or a sphere) 

Cp = heat capacity, J/kg-K (Btu/lb^^-°F) 

h = surface heat-transfer coefficient, W/m^-K (Btu/hr-f t^-°F) 

= heat of fusion, J/kg (Btu/lb^^) 

k = thermal conductivity of solid, W/m-K (Btu/hr-ft-°F) 

s 

r* = radial integration variable, ^j^j^j-gj-face^’^o 
t = ambient temperature, “C (°F) 

t^ = melting point, °C (°F) 

U = overall heat-transfe> coefficient, W/m^-K (Btu/hr-f t^-°F) 

p = density liquid, kg/m^ (lb /ft^) 

pg = density solid, kg/m^ (Ib^^ft^) 

T “ time, hr 
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Two salts were used to compare their lab-scale experimental performance 
with that predicted by the Megerlin solution. The average discharge performances 
of the pure salts LigCOj (mp « 723°C, 1333°F; AH^ “ 6.07 X lo'^ J/kg, 261 Btu/lb) 
and Na^COa (mp « 858°C, 1576°F; AH^ « 2.65 X 10^ J/kg, 114 Btu/lb) were 
selected for this analysis. 

Sol Idlflcatloxt of Pure LipCO^ and NapCO^; Ifovement of the 

Solid-Liquid Interface During System Discharge 

As previously mentioned, T/C's were located at mid-depth in the salt 
at two radial locations [1.8 and 2.3 cm (0.69 and 1.10 in.)] from the HX, and 
their response with time was recorded during discharging. The time required 
for the solidifying front to reach the location of a T/C was determined from 
the traces at the point where the T/G began to indicate a temperature below 
the mp and subsequently yielded a descending cooling curve (characteristic 
of a solid being cooled) . The average times required for the solid front 
to reach the T/C's Isolated from edge effects (such as heat loss) for Lij)C 03 and 
Na 2 C 03 are summarized in Table 19, which shows the radial locations of the T/C*s 
and the calculated Fourier numbers. The theoretical solidification profiles 
determined by numerical evaluation of the Megerlin solution for Li 2 C 03 
(Nph = 0.33, =0.3) and Na 2 C 03 (N^^ = 0.14, Ng^ =0.3) are shown in 

Figures 52 and 53, respectively, with the experimental profile. The observed 
behavior shows the solidification proceeding slower than expected during the 
early stages of discharge and proceeding faster than predicted during the later 
stages, which may be partially attributable to the convective mixing currents 
present in the molten salt. Considering the small size of the unit being 
examined, these numbers agree surprising well with the theoretical model. 


Table 19. COMPARISON OF THEORETICAL AND OBSERVED SOLIDIFICATION 
BEHAVIOR (Li 2 C 03 and Na 2 C 03 ) 


Salt 

System 

Radial 

Location (R-1) 

Average 
Time, min 

K 

Fo 

Observed 

K 

Fo 

Theoretical 

L 12 C 03 

1.76 

17.0 

6.68 

5.8 

Li2C0 3 

3.40 

35.0 

14.00 

19.0 

Na 2 C 03 

1.76 

10.0 

4.07 

3.6 

Na 2 C 03 

3.40 

21.4 

8.72 

10.8 
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The primary source of error for the early-discharge discrepancies is in 
the excess enthalpy of superheat stored as sensible heat in the molten salt at 
the operating temperatures (mp + 50“C). The differences in the later stages of 
discharge are probably a result of the continual radiant heat losses to the ambient. 

It follows from these considerations that the model efficiently predicts 
small as well as larger-scale (engineering units) TES salt performance, as 
found in previous work.^® 

Assessment of Salt TES Applicability 

Because of the complex and diverse trade-offs Involved in selecting a 
salt for a given application and duty cycle, simple comparisons can be 
uncertain, and numbers can often serve only as guidelines when one assesses 
a salt's experimental data for TES applicability. With this in mind, the 
following summaries of salts for solar-thermal applications were prepared, 
based on thermophysical and experimentally determined data. 

• Ll 2 C 03 -CaC 03 55.7-44.3 wt % 

This salt exhibited excellent compatibility with the AISI 304 SS 
cannister and AISI 316 HX tube. Accurate determinations of its heat of 
fusion, thermal conductivity and other thermo physical properties are 
essential, though, for an accurate analysis of its experimental performance. 

It appears to be an attractive PCM for TES applications in the 650° to 700°C 
(1202° to 1292°F) temperature range. 

• NaaCOs-BaCOs 41.8-52.2 wt % 

Reliable thermophysical and transport data are lacking for this 
system also. Excellent long-term (484 hours) compatibility with SS 304 
and 316 containment and stable thermal performance were demonstrated. With 
its higher experimental solidification point of '\/715°C (1320°F) (instead 
of its reported mp of 656°C) , this salt may be utilized in 700° to 750°C 
(1292° t£> 1382°F) applications. 

• (Li-Na-K) 2 C 03 1.21-49.95-48.84 wt % and (Na-K) 2 C 03 50-50 wt % 

The (Na-r.) 2 C 03 binary system exhibited (even during loading) strong 
surface tension creepage properties and was relatively corrosive. Upon 
doping this system with Li 2 C 03 to lower its corrosivity, its heat flux 
increased, though melting somewhat incongruently . Accurate thermal conductivity 
determinations through both solid salts would greatly aid in assessing this 
apparent improvement. 
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• Ll;>CO^ 

Tlu‘ Li;>C 03 unit with tlK> CO;) blanket exhibited the highest heat flux 
of any salts and operated stably for more than 650 hours with no hystersls. 

Good c.ompabllity with the stainless-steel containment was demonstrated, 
although investigations should be conducted into other alloys for long-term 
(30 years) applications. Its high heat of fusion (6.07 X 10‘' J/kg, 261 Btu/lb), 
stable experimental performance, and melting point of 723“C (1339°F) make this 
salt a prime versatile candidate for both large central receiver and small 
dispersed power systems applications in the 700° to 750°C (1292° to 1382°F) 
temperature regime . 

• (Na-K);'C0,, 81.3-18.7 wt % 

This interesting composition requires more in-depth study for one to 
completely understand Its discharge solidification characteristics. Only 
then can the. total impact of its apparent slush behavior on the thermal 
performance and its subsequent applicability be clearly evaluated. As one 
will find in Section 2.'^, though operating 50°C (90°F) higher, generally, 
than the. (Na-K);>C 03 50- ./O wt % system, this salt delivered 31% more, heat 
(Btu/hr-ft^) over the same operating cycle. Its long-term (1032 hours) 

'V ipatibllity with stainless-steel containment and $0. 13/kg ($0. 06/lb) unit cost 
u,,.iin indicate, a strong potential for cost-effective commercial TE8 applicability 
in the broad range of 700° to 800°C (1292° to 1472°F). 

• Na;>C 03 

With its high mp of 858°C (1576°F), medium heat of fusion of 
2,65 X 10*’ J/kg (114 Btu/lb), and extremely low cost of $0. 07/kg ($0. 03/lb), 
this salt seems capable of immediately lending itself to high-temperature 
solar applications in the 800° to 900°C (1472 to 1652°F) range, in which 
the mass of T15S medium is not restricted. Our initial experimental results 
revealed nothing to negate this. They did, nevertheless, clearly illustrate 
the present need to select and test alternate TES containment materials. 

This selection should be performed not only on the basis of salt compatibility, 
but also on the basis of possible interactions with actual working fluids 
(high-temperature air, helium, and molten sodium) or their impurities and on 
the basis of the physical and microstructural stabilities of these alternate 
materials (super alloys) at high temperatures for extended periods (anticipated 
25 to 30-year plant goal lifetimes) . 

• Kg CO 3 

Although it exhibited a high average heat flux of 72,507 W/m*^ (23,000 
Btu/hr-ft^) through two cycles, these numbers are questionable because the 
time at which the HX tube was breached ic unknown. KgC 03 's lower heat of 
fusion (2.0 X 10^ J/kg, 86 Btu/lb) and stec, W decreasing heat flux profile 
result in an uncertain outlook for this salt : ultimate utilization as a TES 
medium in the. 850° to 925°C (1562° to 1697°F) range. Once again, the need 
for investigating new containment materials for these high-temperature 
applications was illustrated. Potassium carbonate should be reinvestigated 
when adequate high-temperature containment materials are identified. 
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• CaCl 2 100 wt % and NaF-KF 32.5-67,5 wt % 

It was found that chlorides and fluorides could not be handled, loaded, 
and operated as easily and safely as carbonates. As previously mentioned, the 
use of either of these salt types on a large scale would require the maintenance 
of a dry, 02 -free inert environment over them to prevent extensive corrosion 
of metallic materials. For the purpose of this investigation, testing of 
chlorides and fluorides, using much the same procedure as that used for 
carbonates, resulted in the loss of containment Integrity and total system 
failure. 


2 . 3 Evaluation of TCE Materials 
Introduction 

An extremely important aspect of TES module performance is the development 
of a means of increasing the average rate of heat transfer between the salt 
and HX, particularly during discharge. Inspection of any salt thermal 
discharge curve reveals a trend of decreasing instantaneous heat flux as the 
thickness of solidified salt between the liquid/solid interface and HX 
Increases. Therefore, the limiting factor for the thermal performance of any 
TES system during discharge is controlled largely by the rate at which the 
latent heat of fusion released at the solid/liquid interface is transported 
to the HX surface, through the growing layer of solid salt. Incorporation 
of a TCE medium within the PCM can significantly influence the performance, 
design, and ultimate cost of the TES system. It allows achievement of relatively 
high average discharge heat fluxes and possibly more uniform thermal input 
to subsequent generators. 

Considerations for the selection of suitable conductivity promoters 
include — 

• Thermal conductivity and heat capacity 

• Compatibility with carbonates 

• Cost 

• Availability 

• Fabrication to desired configuration. 

Experimental Testing 

Table 20 contains a list of potential TCE materials. This table spans 
the range of thermal conductivities of available TCE materials, from 
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21 W/iii-K (12 lUu/hr-n-'’F) for :31(i SS to 308 W/m-K (178 Btu/hr-f t-°F) for 
coppor , 

Of partlculnr noto aro nluinliium (inp 660°C, 1220°F; thermal conductivity 
220 W/m~K, 127 Btu/lir-ft-^F) and copper (mp lOSS^C, 1981°F; thermal 
conductivity 308 W/m-K, 178 Rtu/hr-f t-°F) , both of which poss CHS thermal 
conduc. t IvltlcH 10 tlmea larger than auKtenitlc atalnless steels, and 2 to 6 
times greater than the major constituents of these steels. 

Corros toil 

Alumjm^ Based on the attractive tliermal coiiductivlty 

of aluminum, a compatibility test of a 95 volume percent porous, cmitlnuous 
aluminum matrix (brand named Duocol)* was performed In a (hi-Na-K) 2 CO 3 
ternary eutectic for over 4000 lunirs at 450°C (840°F) . (The mean pore size 
and volume fraction of porosity can be varied over broad ranges during the 
manufacturing of this particular matrix.) 

The cai'bonate melts were contained In alumina crucibles, and a flow of dry 
C0;> gas was maintained over them during testing in Hotpack ovens. Table 21 
summarises the results of this corrosion test. Shown in Figure 54 are an 
unexiKJsed aluminum Duocel cube and one tested for 528 hours. Although one. 
may readily observe the slight discoloraitlon, only negligible weight changes 
occurred after 4180 hours of testing. .Metallography of a polished 
cross-section after testing revealed a very thin, passive corrosion-product 
layer, very likely LIAIO;). 


Table 21. RESULTS OF CORROSION TESTING OF ALUMINUM DUOCEL IN 
TERNARY EUTECTIC 43.5 mol % LioCOj, 31.5 mol % Na;,C0a, 

25.0 mol % K^CO;;,', 450°C 


Sample 
No. 

Exposure 
Time, hr 

... JL. 

Tnitlal Final 

Weight Change 
15 

1 

168 

2.9693 

2.9663 

- 0.0030 

-- 0.1 

2 

528 

3.1034 

3.1014 

- 0.0020 

- 0.06 

3 

1522 

2.9903 

2.9777 

- 0.0126 

- 0.4 

4 

4180 

3.0159 

3.0138 

- 0.002 

- 0.07 


Produced by Energy Research and Generation Inc., Oakland, California. 
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Although it has boon do.inonstrated in IGT’s fviol cell rosoarolr’ that 
aluminum is stable li\ moltoai o,airbonato.s at temperatures remarkably near 
its mp (600°C, 1112°F), aluminum can only be considered as a TCR material 
for low-temperature TKS applications because of its 6bO‘’C (1220°F) mp. 

R eti culat ed Vi tr eous Carbon. Samples of another caiididate 
enhancement material, porous reticulated vitreous carboii (RVC) , wei e also 
chosen for testing. Corrosion tests were, conducted In platinum crucibles 
with the RVC immersed in the 50-50 wt % (Na-K);;C 03 salt at 750°C (1382°F) . 

The initial results indicated imcompatiblllty of the RVC and the salt. This 
was then conclusively determined in a second 400 hour test in which the RVC 
was continuously and totally submerged in the molten salt. The RVC sample 
completely disintegrated, settling into a carbonaceous deposit at the bottom 
of the platinum crucibles. This salt/material incompatibility also eliminated 
RVC from TCE applicability in the high solar-thermal (538° to 871°C) temperature 
range . 

Based on the promising Duocel compatibility results aiid on its high 
thei'inal conductivity, this highly porous body of aluminum was then evaluated 
in a lab-scale TES unit. The results of this test and the evaluation of one 
other TCE concept are presented here, 

Lab- Scale TCE E v aluation 

Aluminum Duocel in (Ll-Na-K) 2 C 03 32-33-35 wt % (L3- 2) . As 
mentioned in Section 2.2, 611 grams (1.3 lb) of the (Li-Na-K) 2 C 03 eutectic x^ere 
initially tested j.n the absence of enhancement materl.al to provide a performance 
baseline to use. when comparing the results obtained with Duocel. Ihis unit 
utilized the original lab-scale design for low temperature, which included 
a slip-fit transite top and tliree salt T/C's shrouded in high-quality alumina 
sheaths. From this initial test, a baseline average heat flux of 6,954 W/m^- 
(2,206 Btu/hr-ft'-) was determined. 

After 10 cycles and 504 hours of operation, a cylindrically fabricated 

block of Duocel x.?as installed into the molten salt, force-fit onto the 1-K 

tube. Figure 55 shows the Duocel after installation. The unit was re-installed 

and the test proceeded in this same cannister an additional 312 hours and 

22 cycles . 

k 

RVC Products Co., Ann Arbor, Michigan. 
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Figure 55. FORCE-FIT INSTALLATION OF DUOCEL 
ENHANCEMENT MATERIAL INTO CANNISTER CONTAINING 
(Li-Na-K)2C03 EUTECTIC 
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Discharge profiles with and without Duocel are compared in Figure 56. 
Although the discharge solidification point remained unchanged ('\^409®C, 

768®P), the graphs are nonetheless significantly different. The time required 
to solidify the ternary out to the T/C at r = 3.02 cm (1.19 in.) decreased 
from 70 minutes to 40 minutes with the porous aluminum Duocel material, a 
43% decrease in solidification time. Also, the unit with Duocel took 33% 
less time to discharge from 450° to 350°C (842° to 662°P) than the baseline 
module, and it reduced the time period of apparent salt temperature thermal 
arrest by 57%. Effective TCE is also illustrated by the fact that Duocel 
caused a 55% decrease in salt radial thermal gradient at the end of discharge 
from 13°C/cm to 29°C/cm (61°F/ln to 133°F/in). This fact is significant 
when one considers that such a I’ncrease occurred across the relatively small 
radius of 1.5 inches. 

The change in the discharge heat flux profiles from no enchancement to 
95 percent porous Duocel is also marked, as shown in Figure 57. The average 
heat flux of 10,135 W/m^ (3,215 Btu/hr-ft^) with Duocel represents a 45% 
Increase in heat conduction through the solid salt to the air working fluid 
over the same operating temperatures. 

After 32 cycles and 816 hours of combined operation, this unit was then 
shut down for post-test evaluation. As shown in Figure 58, the unit was 
sectioned lengthwise, revealing a good salt/TCE-material filling pattern on 
solidification, with little inter-matrix void formation. Visual inspection 
also revealed varying degrees of contact between the Duocel and HX tube, 
varying from direct contact to as much as 2,388 pm (94 mils) separation at 
one point. Despite the force fit achieved during initial installation, the 
Duocel appeared to have relaxed away from the HX during operation. 

One of the primary difficulties in using any reticulated conductivity 
promoter is obtaining a good contact with the HX surface. An alternative 
to force-fitting is to braze the aluminum block to the HX tube. But this 
alternative proved too costly and time-consuming, because further work must 
still be done to develop a braze for use with aluminum and stainless steel 
that is also compatible with carbonate salts. 
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Figure 56. 


(Li-Na-K) 2 C 03 (32-33-35 wt %) DISCHARGE PROFILE 
WITHOUT AND WITH DUOCEL TCE MATERIAL 
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Figure 57. COMPART SON OF (Li-Na-K)2C03 DISCHARGE HEAT FLUX PERFOBMANCE 
WITHOUT AND TOTH CONDUCTIVITY ENHANCEMENT 
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81.3 wt % Na 2 C 03 “ 18.7 wt % K 2 CO 3 (Module Ll-3) . It has been 
determined from previous TES work performed at IGT^® that heat transfer 
resistance in the PCM is the major TES performance-limiting factor. This 
effect, which results from the low thermal conductivity in the solid phase, 
reduces the heat flux through the solid PCM that forms on the salt-side 
heat-transfer surface. It is the primary limitation to using a passive 
HX in the TES system. With these problems in mind, this particular (Na-K) 2 C 03 
system was chosen as the third TCE concept to evaluate. From tnis system's phase 
diagram, a mixture of 85-15 mol % Na 2 C 03 -K 2 C 03 appeared to represent a system 
that would solidify incongruently over the range 790° to 740°C (1454° to 
1364°F) upon discharging. 

It was considered that such incongruent solidification could result 
in a "slush- type" salt behavior, perhaps eliminating the formation of a 
continuous crust of solid salt on the HX surface. In this way, heat transfer 
could be enhanced by convection in the remaining dispersed liquid phase. 

To this end, 821 grams (1.8 lb) of this salt miocture was loaded and run 
as both a prospective high-temperature solar-thermal TES medium and TCE 
material. Examination of the salt temperature profile as a function of 
time (Figure 59) reveals that no actual thermal arrest occurs in the temperature 
range 790 to 740°C (1454° to 1364°F) ; instead, two slight inflection points are 
exhib’‘Ued in this area. This apparent behavior is consistent with the incongruent 
character of the published phase diagram. 

Also, from Run 21 (Figure 60), the average heat flux of 48,375 W/m^ 

(15,345 Btu/hr-ft^) was calculated. This represents a 31% increase in heat 
flux to the air working fluid by this mixture over that of Module Ll-1, 
containing 50-50 wt % Na 2 C 03 ~K 3 C 03 over the same "working cycle." Also, 
throughout its experimental life, this system's performance was consistent 
and stable. 

This salt system was subjected to prolonged cycling of 38 charge/discharge 
runs and 1032 hours at operating temperatures of 813 to 713°C (1495° to 1315°F), 
with no system failures or ‘^alt leakages. After a forced cooldown at 5.1 m^/hr 
(180 ft^/hr) , the cannister was sectioned at mid-salt. 


117 


INSTITUTE 


0 F 


GAS 


technology 


SALT TEMPERATURE VS. TIME 




Oo ‘3dnivd3dW3i i3nino aiv 



0 « ‘3dfUVd3dV^31 nvs 


GAS 


INSTITUTE 


0 F 


TECHNOLOGY 


'igura 59. TYPICAL DISCHARGE PEREOEMANCE OF (Na-K> 2 C 03 SLUSH 
AT 5.1 m3 /hr (180 ft3/hr) AIR FLOW RATE 
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Figure 60. LIFETIME DISCHARGE HEAT FLDX OF (Na-K)^CO- SLUSH 
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Some interesting observations come to light when this salt’s crystal 
patterns (Figure 61) are compared to the solidified surface of the similarly 
cooled Na 2 C 03 -BaC 03 congruently melting eutectic of Module L3-3. Rather than 
the pattern of crystals Indicative of rapid cooling (Figure 47), one encounters 
unoriented crystals of similar size and geometry occurring from the HX tube 
to the can inside-diameter, 

In view of these facts, the thermal conductivity enhancing "slush" 
behavior of the (Na-K) 2 C 03 mixture is notably confirmed by the — 

• Absence of thermal arrest in its salt temperature discharge profile 

• Incongruent solidification patterns 

• 31% Increase in average heat flux over the 50-50 wt % congruently 
melting (Na-K) 2 C 03 eutectic. 

The post-test chemical analysis of this salt, which revealed no decomposi- 
tion after 1032 hours and little corrosion, is summarized in Table 22. 

Table 22. (Na-K) 2 C 03 81.3-18.7 wt % "SLUSH" PRE- AND POST-TEST CHEMICAL 

ANALYSIS 


Chemical Analysis, Wt 

% 

Na 

K 

CO 3 

Fe 

Cr 

Ni 

Pre-Test 


35.27 

10.58 

54.16 

<0.01 

<0.01 

<0.01 

Post-Test (r = 0.23 

in. ) 

35.16 

10.02 

54.59 

0.05 

<0.01 

<0.01 

(r = 0.68 

in . ) 

35,00 

11.65 

55.35 

0.02 

<0.01 

<0.01 

(r = 0.95 

in . ) 

35.07 

11.99 

55.15 

0.01 

<0.01 

<0.01 


Techno-Economic Consideratlona of TCE Materials 

When considering the feasibility of TCE conceptSj. one must evaluate both 
the technical and economic aspects of the problem, that is, the cost/performance 
trade-offs. In costing an actual TES Stt'b8y8tem)i the major component costs to 
be considered include “ 

• Salt PCM 

• Containment and HX materials 

• Fabrication and assembly 
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• Insulation 

• Instrumentation and controls 

• TCE materials or active HX system. 

Economic analyses indicate that TES subsystem costs are dominated by 
salt and containment/HX costs. Yet the overall subsystem cost may be either 
salt-cost intensive or HX-cost intensive, depending on the application's duty 
requiremenls, which dictate selection of salt properties and HX configuration. 
Although TCE materials do indeed add an incremental cost, such an approach 
may reduce the required heat transfer area in HX-cost-intenslve cases enough 
to result in a net decrease in total TES subsystem cost by reducing HX 
materials and fabrication costs. Other items to consider when one is selecting 
TCE materials and designing subsystems are the sensible heat storage capability 
added to the subsystem and the nominal increase in required TES volume caused 
by displacement of salt and its latent heat. 

Conductivity Enchancement Factor 

The most practical porous aluminum densities for TES application are in 
the range of 5% to 20% of theoretical density (i.e., 0.95 to 0.80 volume 
fraction porosity). In this range, the porosity is essentially completely 
interconnected, so that both the salt and aluminum enhancement form continuous 
phases throughout each segment of the TES. To a first approximation, then, 
the effective thermal conducitivity of the two-phase mixture can be expressed 
as — 


ef f 


f k + f.-k.- 
s s A1 A1 


- 


* ^aAi 


where ^ 


f = volume fraction of aluminum enhancement ;material 
A1 

f = volume fraction of salt 
s 

^A1 ~ thermal conductivity of aluminum enhancement material 
at theoretical density 


kg = thermal conductivity of solid salt. 


(32) 
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Till' inamirai'turi'rs o! l:lu> Ouocol inatorlal expor linoutally moasui'i'd tho Llu'rmnl. 
I'lnuiuct ivl Ly of poroua aluminum as a function of donsity and obtained the 
solid eiu've shown in Figure b2. This curve can be compared Vi'lth the calculated 
dashed curve, which Is based on Kiiuatlon 12 and neglects any contribution to 
conductivity by the air in the pores, in tlie high-porosity range investigated, 
tho term overpredicts the actual conductivity of tlie porous aluminum. 

Oorrectlng for the difference between measured and calculated values, 
liquation can be modified to 


I'ff 


( 1 - 


f.,')k + 0.5(1 

A1 H 




A tliermal conductivity enhancemeiU. factor, F, was defined as the ratio of 
effective conductivity of the solid salt/porous aluminum composite material, 
to till' conductivity of solid salt alone, as sliown in liquation 3A; 


k 


off 

k 

s 


( 1 - 


•f. ,)k 
A1 s 


0.5(1 


Cl^) 


liquation was used to calculate values of for a series of assumed 

values of k^ and f^^j , and tlie results appear in Figure (i3. For salts liaving 
k values in the range of 0.87 to 8.(i5 W/m-K (0.5 to 5 lUu/hr-ft-'’F) , use 
of a high porosity aluminum enhancement material would be expected to result 
in an increase In k^^j.|., corresponding to approximately 1.12 W/m-K (0.(i5 Btu/ 
lir-ft-”]'’) for each percent of theoretical density of tlie aluminum used. For 
a given aluminum density, however, tlie largest degree of enliance.ment 
(larger F) is realized for salts of lowest conductivity. 

Cost , Performaiu'.o and Design. T.nqilicatluns 

Use of a TCli material will impact on an actual TK8 system in several 
areas that must be evaluated: 

• Thc.rease in the effective thermal conductivity, through 

the solid salt layer 

• Decrease in the recgiired liX surface area 

• Supplement of sensible heat storage cnpacity to the TK8 system 

• Increase in required system volume, a result of displacement 

of the. salt 

• Increase in capital costs associated \'/lth TOli medium. 
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Additional costs are undeniably associated with any enhancement material. 

It must be realized that certain counterbalancing trade-offs do exist: the 
increased will also lead to a need for less HX area. This may significantly 

offset the cost of enhancement when the net TCE cost is viewed as being 
approximately equal to the cost of enhancement material minus the cost of 
savings of reduced HX area, considering the salt storage volume displaced 
by TCE is negligible. 


As a first-cut estimate of the performance-cost trade-offs involved in 
the use of conductivity enhancement materials, it is useful to consider the 
relationship between the enhancement factor, F (F = 

capital cost of the materials. F factors based upon the curves of Figure 63 
have been calculated and are plotted in Figure 64 as the ratio of enhancement 
factor to unit price, F:$/m^ (F:$/ft^), as a function of percent 
theoretical density of aluminum for a range of assumed values. These 
curves are based on a theoretical density of aluminum of 2,712 kg/m^ (169.3 
Ib/ft^) and an estimated bulk-*quantity cost for the fabricated porous 
aluminum struture of twice the raw material cost [i,e., 2 X $l*l/kg ($0. 50/lb); 
see Figure 65 for assumed cost versus density]. As expected, the enhancement 
factor-to-cost ratio decreases with increasing k^ values. It also decreases 
somewhat with increasing aluminum density, although this trend mostly reflects 
the assumed shape of the cost versus density curve, shown in Figure 65. 


The effect of an enhancement material on performance is best shown 
through use of the heat transfer model based on Megerlin’s solution to the 
moving boundary problem. 

Parametric curves previously derived for the engineering-scale units, 

which contain LiKG 03 and are cooled by high-velocity ambient air, are shown 

in Figure 66. These curves show dimensionless time required to solidify 

the LiKC 03 salt of Engineering Unit 2 to a dimensionless thickness (R-1) . 

For these units, - 0.57 and = 0.76, in the absence of any conductivity 

enhancement. Approximately 5 hours are required to discharge these units; 

this time corresponds to = 20.7 before correcting for removal of superheat 

prior to the start of solidification. Such a discharge results in a solid 

salt layer of reduced thickness, R-1 - 3.7, and appears as Case IV in Figure 

66 for k - 2.25 W/m-K (1.3 Btu/hr-f t-“F) . But if porous aluminum 
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Figure 64. ENHANCEMENT FACTOR-TO-COST RATIO VERSUS DENSITY OF POROUS ALUMINUM 

FOR SEVERAL ASSUMED VALUES OF k 
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enhancement of 8 % density were introducted into the system, would 

be expected (by Equation 33) to increase from 2.25 W/m-K (1.3 Btu/hr-f t-°P) 
(k for solid LiKCO? without enhancement) to 11.2 W/m-K (6.5 Btu/hr-f t-°F) , 
and an F-factor of 5 (Case VI) . Figure 66 also shows that the parameters N 
and Np^ have different dependences on 

1 


Bi 


^Bi “ k 


off 


and 


N k 

eff 


so that now - 0.15 and 103. As indicated, this 5-^hour discharge 

with an enhancement factor of 5 would allow solidification of a layer of 
thickness R~1 = 5.3. With the relative solidified volume being proportional 
to r^ (from R = r/r , r - 1 in.), the enhancement material could increase 
the average heat flux so that approximately twice (1.8 X) the salt volume 
could be solidified in the same time. In a large TES unit with multiple 
ex’s, this is equivalent to reducing the required HX surface area by a factor 
of two, and it results in corresponding cost reductions. 

Because of this potential for increasing the thermal performance and/or 
decreasing the cost of TES substyems by incorporating TCE materials, it becomes 
apparent that an extensive search for and evaluation of such materials for 
use in molten carbonate environments should be extended to the 700° to 870°C 
(1300° to 1600°F) solar-thermal regime. 

2.4 Conceptual Design of TES Module for Solar-Thermal Power System 

Use of the Heat Transfer Model in Large-Scale TES Solar-Thermal 
Power System 

Introduction 

The phase-change heat-transfer model developed by Megerlin and modified 
during our previous TES program at IGT can also be used in large-scale system 
design. The existing model can be used to determine the thickness of 
solidified salt around a cylindrical internal HX tube that carries a working 
fluid as a function of the time it has cooled from its mp. This 
solidification is measured in terms of the radius of the internal well. 
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For a conceptual system design that has a large molten salt bath with 
a number of pipes parallel to each other in a well-formed matrix, carrying 
the working fluid through the salt bath, the model will determine the pipe 
spacing, the number of pipes required, and the length needed (given the 
pipe radius) for a particular storage requirement. Thus, the energy discharged 
from a unit is related to the physical design of the system, allowing a 
balance between the energy required by the working fluid and the energy that 
is discharged as a result of the system design. 

TES Sub-System for Large-Scale Solar-Thermal Power PI nt (Air Fluid ) 

The heat transfer model has been applied to the scale-up design of a 
storage sub-system of a central-receiver solar-thermal power-plant that 
utilizes a high-temperature Bray ton cycle with high-temperature, high-pressure 
air as the heat transport fluid. 

Such a central-receiver solar power-plant concept has been proposed by 
Boeing,^® utilizing a high- temperature, closed, Brayton-cycle thermal engine, 
with helium or air as the heat transport fluid. Figure 67 is an overall view 
of the central receiver and hellostat field, with the power plant and storage 
subsystem at the base of the receiver tower. 

Figure 68 illustrates a solar-thermal power plant schematic incorporating 
the TES unit. The primary flow circuit from the receiver is through the 
turbine(s), the recuperator (low-pressure side), the precooler, the compressor (s) , 
the recuperator (high-pressure side), and back to the receiver. The alternative 
thermal-energy stroage loop extends solar-plant electrical-power production 
into non-insolation hours. The resultant phase-change TES conceptual design 
is summarized in Figure 69. The phase-change medium is contained in an 
Insulated underground vessel. The working fluid is channeled through the 
specified design number of straight superalloy tubes and collected by manifolds 
at the top and bottom of the vertical vessel. The carbonate salt can then be 
contained as a bath immersing the HX tubes. The vertical flow arrangement, with 
the hot working fluid entering the top manifold during charging and the cold 
fluid entering the bottom manifold during discharging, will maintain the solidifed 
salt below the liquified portion. This arrangement will provide for unrestrained 
volume expansion during salt melting, reducing stresses on the system 
containment and HX tubes. 
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Figure 69. PHASE-CHANGE TES CONCEPT 
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The system is designed to deliver 50 over a 6-hour period, thus 

requiring a system capacity of 1.08 X 10^^ J (1.02 X 10^ Btu) . 

Conceptual design of the TES subsystem is based on use of pure Li 2 C 03 
(mp of 723°C, 1333“?) as the PCM for the purpose of heating the 3.45 MPa 
(500 psia) air from 538°C (1000°F) at the inlet to 704°C (1300“F) at the 
outlet. 

For a 50 TES delivery, the air flow rate then must be — 

m^ = 50 MW X (1.08 X 10 J/MW) ( (1/1.87 X 10^ J/kg) = 

= 9.66 X 10^ kg/hr = 2.13 X lO® lb /hr (35) 

m 

Also, from the thermal capacity — 

m = (1.08 X 10^2 j) (1/6.07 X 10^ j/kg) 
s 

= 1.78 X 10® kg - 3.93 X 10® lb (36) 

“ m 

of salt (heat of fusion only) are required for storage. The volume of salt 
that must then be solidified (using the room temperature density for a first 
approximation) is — 

= (1.78 X 10® kg)/(2108 kg/m3) = 843.8 m®= 2.98 X 10^ ft^ (37) 

This now determines the "fixed volume" requirement that the heat transfer 
model will be applied to solve. The phase-change number determined for a 
mean air temperature of 621°C (llbO^F) is — 


Nph 


^p(s) ^*'m *"a^ ^s 


= 1.94 = 2 


(38) 
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where 

« 6.07 X lO''^ j/i^g (261 Btu/lbJ 

» 1834.1 Icg/m^ (114.5 Ib/ft^) 

p = 2108 kg/m^ (131.6 Ib/ft^) 
s 

= 2625.1 J/kg-K (0.627 Btu/lb-°P) 

= melting point temperature = 723°C (1333°!') 
t = mean air temperature = 621. 1“C (1150°F) 

cl 


The Megerlin plots for N 
Figure 70. 


ph 


= 2 and a range of Biot numbers appear in 


In this example, pipes 3.81 cm (1.5 in.) in diameter were chosen; the 
number of pipes required, their spacing, and their length had to be determined. 
The Fourier number corresponding to a 1,9 cm (0.75 in.) radius and 6-hour 
discharge time (t) is — 


N. 


Fo 


a T 

^ 2 
O^ 


= 15.67 


where a = 


k 
c p 

P s 


(39) 


Given the pipe diameter, (d = 3.81 cm, 1.5 in.) 2r^, the Biot number is 
specified for using the Megerlin plots in Figure 70 to determine r^ , which 
will yield the thickness of solidified salt (rj - r^) and the pipe spacing 
(2ri). The pipe length is determined from the fixed volume requirement — 


V 


Li = 


1 - (r 2- 
1 


r 2) n'lr 

o 


43.8 m^/(ri^- r ^) nrr 


= [2.98 X 10‘^ ft^/(rj2 - r^^) mr] 
With the Biot number being given by — 


N„_. 


hr 
f o 


Bi k 


(40) 


(41) 
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an expression for the air-side heat transfer coefficient, h^, must be derived. 
In general, one may assume fully turbulent flow in the air tubes. This 
assumption thus permits use of the Nusselt equation ~ 


= 0.023 

Nu Re 


N 


0.2 

Pr 



(42) 


to derive the expression for h^. Solving Equation 42 for hf yields ~ 


hf “ 


0 ? 

0.02 kr 

Pr f 


m. 


•1 0.8 




> —1 • 8v y -0 . 8v 

(r^ ) (n ) 


(43) 


where n = number of pipes with radius r . Substituting into Equation 43 the 
appropriate properties for air at 621°C (1150°F) and 3.45 MPA (500 psia) , while 
setting r^ = 1.9 cm (0.75 in.) and = 9.66 X 10^ kg/hr (2.13 X 10® Ib^^/hr) 
yields 

h^ = 47,363/n®’^(W/m-K) [27,366/n^*® (Btu/hr-ft -°F)] (44) 


It can be seen from Equation 44 that and r^^ are both dependent on 
the number of pipes to be used. The pipe length, can be determined for 
any number of pipes by selecting n and calculating h^ from Equation 44, 
determining and using it to find r^^ from the Megerlin plots, then 
substituting r^ and n into Equation 40, and solving for 

Wlien one has calculated a pipe length consistent with a given number of 
tubes, n, that are sufficient to remove the heat from the salt volume, he 
must consider if sufficient heat-transfer surface is present to heat the air 
to the required outlet temperature. 

Assuming that the PCM remains at its melting point, the enthalpy change 
in the air working fluid as it progresses through the pipe is described by — 


^^Vd S " '"f ^*"m “ t) = hf (726‘*C - t) [hf (1338. S^F - t)] (45) 
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Rearranging Equation 45, integrating, and evaluating [given that at L = 0, 
t ^ 538“C (1000°F), and at L » L2. t “ 704°C (1300'’F)], leads to - 


^ n IT he d 

f 


(46) 


The total heat transfer requirement will be satisfied when n is chosen so 

that Lj = Lg* thus meeting both the fixed-volume and air-temperature requirements. 

The solution is obtained graphically by using Equations 40, 44, and 
46, and the Megerlin plots of Figure 70 to derive Table 23. 

The values of Li and L2 are shown as a function of n in Figure 71. The 
intersection of the two curves, representing the design solution, occurs at 
n = 52,000 and Lj = L2 = 3.38 m (11.1 ft). This design operates with hf = 

7.96 W/m-K (4.6 Btu/hr-f t^-"F) and r^ = 4.34 cm (1.71 in.). Therefore, the 
design conditions for the TES operation can be met by solidifying 1.8 X 10® kg 
(3.9 million lb) of 112003 cylindrically around 52,000 pipes, 3.81 cm (1.5 in.) 
in diameter and 3.38 m (11.1 ft) long, spaced 8,7 cm (3.42 in.) apart, 
over a 6-hour period. A top view of this TES subsystem, with hexagonally 
arrayed tubes, appears in Figure 72. As a final check, the total area open 
to air flow is — 

A = mr vj = (52,000) (ir) (f^)^ = 59.3 m^ 

^ (52,000) (tt) (^^)^' = 638.1 ft^ (47) 


The linear air velocity is -- 


y (9.66 X 10^ kg/hr) (hr/3600 sec) (0.0699 m^/kg) « q 32 m/sec 

59. 3 

_ (2.13 X 10® Ib/hr) (hr/3600 sec) (1.12 ft®/lb) , 

_ 638.1 ft^ ■ r»^/sec 


(48) 


Where the specific volume of air at 3.45 MPa (500 psia) , 621“C (IISO^F) is 
0.0699 m®/kg (1.12 ft®/lb), the Reynolds number is ~ 


/ (1.5/12 ft) (1.04 ft/sec) (0.891 lb/ft3)(3600 sec/hr) 

"Ee ■ ■‘''Pf/lJf" ■ 0.0466 lb/£t-hr 


= 4.65 X 10®/n = 8,945 (n = 52,000) 


(49) 
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which confirms the turbulent flow assumption and verifies the use of the 
Nusselt equation to derive the air-side heat transfer coefficient. 

It is to be kept in mind that while the above design does Indeed satisfy 
the performance requirements, it is not necessarily the optimum design. An 
optimization analysis must Include studies of the effects of pipe sizing and 
of the use of conduction enhancement and determination of the cost-effectiveness 
of each design, which is beyond the scope of this project. Nevertheless this 
example does show the importance of the heat-transfer model and does identify 
the key parameters involved in designing a latent-heat TES subsystem for 
solar-thermal applications. 
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CONCLUSIONS 

Steam-Electric Power System Deflrtitlon 
Xu this project, we. — 

• Defined a TES utilization concept of an electrical generation subsystem, 
composed of a multistage condensing steam turbine and a TES subsystem 
with a separate povfer conversion loop. [The primary plant provides 
steam at 538°C (1000®F) and 16.5 MPa (2400 psla).] 

• Evaluated conceptual designs for a 100 TES peiiking system providing 

steam at 316°C (600"F), 427°G (800°F) , and 454°C (SSO^F) at 3.79 MPa 
(550 psia). 

316 ° C case: Three stages; charge-to-discharge mass ratio of 

1.12 with a 20% system efficiency. Average break-even storage 
cost, with electrical generation subsystem cost of $200/kW 
and desired energy costs of 60 mills/kWhr is $43.5/kWhr^^^ . 

42 7 °G ca se: Six isothermal sect. ions; charge-to-discharge 

mass ratio of 1.2 with a 28% system efficiency. 

45 4°C cas e; Seven isothermal stages; charge-to-discharge. 
mass ratio of 1.3 with a system efficiency of 30%. Average 
break-even storage cost is $52/kWlir^^j . 

• Goncluded TES-peaklng-system utility applications that employ a separate 
power-conversion loop appear technically viable when interfaced with 
high-temperature, high-pressure base-load steam-electric plants. 

« Concluded the economic, feasibility of the TES peaking system depends on 
the development of reasonably priced containment vessels/HX's with 
high heat flux rates. 


Salts 


• Carbonate salts are attractive candidates for thermal energy storage 
applications . 

-- Tliey have excellent thermal, chemical and physical properties 
atid stabilities. 

— Safe, simple handling procedures can be used. 

• The feasibility of using alkali/alkaline earth carbonates as latent-heat 
thermal-storage materials over the temperature range 397° to 898°C 
(747° to 1648°F) was demonstrated. 

• L 1 KC 03 (mp 505°C, 94l°F) exhibited excellent chemical, physical 
and thermal stability for over 5600 hours and 130 cycles and was 
compatible with Type 316 SS. No change in salt chemistry occurred. 
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• Pure U 2 CO 3 [6.07 X 10® J/kg (261 Btvi/lb), $2. 05/Kg ($0. 93/lb)] and Na 2 C 03 
[2.63 X 10^ J/kg (113 Btu/lb), $0. 07/Kg ($0. 03/lb)] appear attractive 

as PCM’s in approximate temperature regimes of 700“C (1300®F) and 815“C 
(1500“P), respectively. 

• The Na 2 C 03 -BaC 03 mixture is an attractive, low~cost ($. 18/Kg, $0.08/lb) 
mixture for ^i700®C ('vl300“F) applications. 

• Salt of 85 mol % Na 2 C 03 - 15 mol % K 2 CO 3 solidified incongruently over 
the temperature range 785® to 740®G (1445® to 1364®F) . 

• Salts of CaC 03 and BaGOs were chemically stable as mixtures with alkali 
carbonates over the 650® to 700°G (1200® to 1300®F) range. 

• Addition of 15 mol % GaC 03 to (Li-Na“K) 2 C 03 ternary eutectic lowered 
the melting point from 397° to 383®C (747® to 721°F) . 

• Garbonates can be melted and loaded in air environment safely and 
acceptably. 

• Chlorides and fluorides require more sophisticated handling fabrication 
and operation procedures. 

- Vapors are toxic in air. 

— They must be handled and cycled in dry, inert environments to avoid 
accelerated containment corrosion by oxygen/raoisture contamination. 

Containment and HX Materials 


• Types 304 and 316 SS demonstrated good compatibility over a range of 
370° to 650°C (700® to 1200®F) . Salt-side oxidation of Type 316 SS 25 pm 
(1 mil) after 5650 hours in contact with LIKCO 3 at 482° to 524°C 

(900° to 975°F) occurred. 

• Prolonged operation above 650°C (1202°F) results in greater salt creepage 
and increased oxidation, carburization, grain boundary, carbide precipi- 
tation and brittle sigma-phase formation in austenitic stainless steels. 
Use of superalloys sould be investigated. 

• Use of CO 2 blanket over molten carbonates above 650°C reduces the 
rate of oxidation and carburization of stainless-steel materials. 

TCE Materials 


0 Porous aluminum material is attractive in carbonates to near its mp of 
^650°C. It demonstrated long term stability at 450°C (840°F) . 

0 Addition of 5 vol % aluminum Duocel increased the average discharge heat 
flux in ternary (Li-Na-K> 2 C 03 eutectic (mp 397®C, 747®F) by 'V'45%. 
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• Reticulated carbon matrix is not compatibile with carbonates under air 
at 750°C (1382°F). 

• An iiTcmigruent slush concept of 85 mol % NaoCOs - 15 mol % R 2 CO 3 salt 
exhibited a heat flux 31% higher than the 45-56 mol % (Na-K) 2 C 03 system, 

• Cost/benefit trade-off of TCE materials for solar-thermal systems depends 
on system design and duty cycle, and salt and TCE material characteristics. 
Further work is required for applications above 650°C (1202®F) . 

TES Module Thermal Performance 

• Over 1000 hours of stable laboratory -scale and 5000 hours of stable 
engineering-scale performance have been demonstrated with carbonates. 

• Repeated thermal cycling (20-130 cycles) did not affect physiochemical 
properties or cyclic thermal performance of any carbonate salts studied. 

• Baseline thermal discharge characteristics and heat fluxes for salts 
considered for solar-thennal applications [37,830 - 75,660 W/m^ (12,000 - 
24,000 Btu/hr-ft^) in 0.1 KWlirth units] were established. 

0 The heat-transfer model based on Megerlin's solution shox^s agreement 
between predicted and observed interfacial location during heat 
removal from PCM. 

• The heat-transfer model for scale-up to multi-tube 50 HX 

configuratioxi at 704°C (1300°F) was utilized. 

• The thermal enhancement concept provides more efficent utilization 
of heat-transfer area. 
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RECOMMENDATIONS FOR FUTURE RESEARCH 

Screen and evaluate high-temperature alloys as carbonate containment 
materials for 650° to 870°C (1200° to 1600°F) applications. 

Pursue development of cost-effective salt, containment, and insulation 
materials . 

Scale up to 10 KWlirth modules and conduct endurance testing of the 
most promising salts identified in lab-scale (0.1 kWlirj.jj) studies, 
utilizing working fluids considered for commerical applications. 

Continue cost/benefit trade-off analyses and compatlbility/performance 
studies on TCE materials in a 650° to 870°C (1200° to 1600°F) temperature 
range. 

Further develop and experimentally verify the heat-transfer model to 
account for TCE materials, high- temperature working fluids, alternate 
HX configuration, and superheating and convective mixing effects on 
solidification and melting behavior, 

Maintain cognizance of and active interaction with ongoing TES 
systems design and analysis studies to provide effective feedback 
to salt/HX R&D efforts. 
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Non-SI Unit 

Table A-1. CONVERSION FACTORS 
Operation 

SI Unit 

atm 

X 1.013 X 10^ 

Pa 

Btu 

X 1055.06 

J 

Btu/lb 

X 2326.00 

J/kg 

Btu/lb-°F 

X 4186.80 

J/kg-K 

Btu/hr-ft-°F 

X 1.731 

W/m-K 

Btu/hr-ft^-°F 

X 5.678 

W/m2-K 

Btu/hr-ft^ 

X 3.154 

W/m^ 

Cal 

X 4.18 

J 

°C 

+ 273.15 

K 

cP 

X 0.001 

Pa-s 

“F 

(5/9) (°F-32) 

°C 

ft 

X 0.3048 

m 

CM 

X 0.0929 

m^ 

ft2/hr 

X 0.0000258 

m^/s 

ft3 

X 0.0283 

m3 

in. 

X 0.0254 

m 

in. 2 

X 0.00064516 

m^ 

lb 

X 0.4536 

kg 

Ib/ft^ 

X 16.018 

kg/m3 

Ibm/ft^-s 

X 4.882 

kg/m^-s 

mil 

X 0.0000254 

m 

psia 

X 6894.8 

Pa 
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APPENDIX B. Carbonate Phase Diagrams 
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Figure B-3. PHASE DIAGRAM FOR T1 
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